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ABSTRACT 


Previous phylogenetic analyses of morphological and rRNA data indicated that Gnetales are the closest living 
relatives of angiosperms but gave different basal angiosperm relationships. A two-step morphological analysis of seed 
plants (including fossils) and angiosperms rooted the latter near Magnoliales, with tricolpate eudicots and paleoherbs 
(herbaceous magnoliids and monocots) forming a clade, whereas analyses of rRNA sequences rooted angiosperms 
among paleoherbs, with eudicots and woody magnoliids forming a clade. Experiments with a revised seed plant 
morphological data set raise further questions: when angiosperms are scored like different angiosperm subgroups, 
they associate with different outgroups, although Gnetales are their closest living relatives. To test whether morphological 
and rRNA data are seriously contradictory or rather complementary, with inconsistencies being a function of better 
resolution in different parts of the tree, we experimented with morphological and rRNA data sets including the same 
six extant "gymnosperm" and 12 angiosperm taxa. Both analyses again associate angiosperms and Gnetales. The 
morphological analysis differs from previous ones in placing Nymphaeales and monocots at the base of the angiosperms, 
but trees rooted next to Magnoliales are only one step less parsimonious. As in previous studies, the rRNA analysis 
roots angiosperms next to Nymphaeales and breaks up the eudicots. Bootstrap and decay analyses of the rRNA data 
show strong support for the monophyly of angiosperms and Gnetales and their sister group relationship, but low 
support for groupings within angiosperms. However, one or another group of paleoherbs is basal in most bootstrap 
trees, A combined analysis favors a paleoherb rooting, but other relationships agree with the morphological results; 
particular, eudicots form a clade. The conclusion that Gnetales are the closest living relatives of angiosperms 
permits a wide range of morphological scenarios, depending on the arrangement of fossil outgroups. Discovery of 
fossils on the long branch leading to angiosperms, methods of factoring out artifacts in rooting, or molecular data on 
the control of floral morphogenesis in angiosperms and Gnetales may be required for further progress in unraveling 

origin of angiosperms, 


БЕРЕ SSS 


Р Recent cladistic analyses of morphological and Our initial morphological analysis of seed plants 
ribosomal RNA data have led to apparently con- (Doyle & Donoghue, 1986), which included both 
flicting results on the position of angiosperms among living and fossil taxa, treated angiosperms as а 
seed plants, basal relationships among angiosperms, single taxon, modeled on Magnoliales and Winter- 
and resulting scenarios for the origin of angio- асеае, following the consensus at that time. This 
‘perms (Crane, 1985; Doyle & Donoghue, 1986, study was designed especially to test the previous 
1992; Donoghue & Doyle, 1989a; Zimmer et al, analysis of Crane (1985) and other current ideas 
1989; Loconte & Stevenson, 1990, 1991; Taylor оп seed plant phylogeny, by including both char- 

Hickey, 1992; Hamby & Zimmer, 1992). The acters used by Crane and conflicting similarities 
Purpose of this paper is to investigate the causes that he had omitted. Contrary to our expectations, 
and significance of these conflicts, and what if any but as Crane had found, this analysis indicated that 
robust conclusions concerning the origin of angio- seed plants are a monophyletic group, with conif- 
‘perms can be drawn from these data. eropsids (cordaites, conifers, ginkgos) derived from 
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advanced “‘seed ferns” with platyspermic seeds and 
saccate pollen (roughly as proposed by Rothwell, 
1982), and that angiosperms belong to an “ап- 
thophyte" clade also including Bennettitales, Pent- 
oxylon, and Gnetales, which is in turn nested among 
so-called Mesozoic seed ferns (corystosperms, glos- 
sopterids, Caytonia). In response to a cladistic 
analysis involving only extant seed plants by Lo- 
conte & Stevenson (1990), we revised this analysis 
with generally similar results (Doyle & Donoghue, 
1992), except for more definite placement of cy- 
cads among platysperms and more uncertainty on 
the position of ginkgos, which may be associated 
with the Permian-Triassic seed fern Peltasper- 
mum instead of conifers (cf. Meyen, 1984). 

A subsequent analysis of extant angiosperms 
(Donoghue & Doyle, 1989a), using the results of 
the seed plant study to polarize characters within 
the group, yielded trees rooted in or next to the 
"woody magnoliid” order Magnoliales (specifically 
families with granular exine structure, which ex- 
cludes Winteraceae and Austrobaileya). This re- 
sult is generally consistent with conventional views 
on angiosperm evolution (e.g., Takhtajan, 1969; 
Cronquist, 1981, 1988). The remaining angio- 
sperms form four major clades: two other groups 
of woody magnoliids, (1) Laurales, including Chlo- 
ranthaceae, widely discussed because of their un- 
usually simple flowers, and (2) winteroids, including 
Winteraceae, Illiciales, and possibly Canellaceae; 
(3) dicots with tricolpate and derived pollen, later 
designated eudicots (Doyle & Hotton, 1991); and 
(4) "paleoherbs," consisting of herbaceous or semi- 
herbaceous magnoliids (Aristolochiaceae, Lactoris, 
Piperales, Nymphaeales) and monocots. Somewhat 
similar results were obtained by Loconte & Ste- 
venson (1991); the most important difference, 
placement of Calycanthaceae and Idiospermata- 
ceae at the base of the angiosperms, may be due 
to their use of Recent plants only as outgroups, 
since the closest outgroup, Gnetales, shares fea- 
tures such as opposite leaves and two-trace nodes 
with Calycanthales. 

In contrast, analyses of partial 18S and 26S 
rRNA sequences from angiosperms and other ex- 
tant seed plants (Hamby & Zimmer, 1992), chosen 
to provide a test of current morphological hypoth- 
eses, placed a series of paleoherb taxa at the base 
of the angiosperms (Nymphaeales, Piperales, and 
monocots). Woody magnoliids and eudicots formed 
a more derived clade, within which detailed rela- 
tionships are poorly resolved. These results give a 
different picture of the first angiosperms: they would 
be herbaceous or nearly so, with palmately veined 
leaves and anomoc ytic stomata, rather than woody 


with pinnately veined leaves and paracytic stomata. 
This recalls the views of Burger (1977, 1981) on 
the primitive status of Piperales and/or monocots, 
although not in detail. 

More recently, the view that the first angio- 
sperms were “‘paleoherbs’’ was elaborated by Tay- 
lor & Hickey (1990, 1992), based on recognition 
of an Early Cretaceous (Aptian) paleoherb fossil 
from Australia and their own morphological anal- 
ysis. Their tree differs in that Chloranthaceae are 
basal, followed by Piperaceae, which contrasts with 
the position of Chloranthaceae in Laurales in the 
analysis of Donoghue & Doyle (1989аг), well re- 
moved from Piperales, and in the woody magnoliid- 
eudicot clade in Hamby & Zimmer (1992). This 
difference may be partly a result of treating cor- 
related ovule features related to orthotropy (which 
was assumed to be primitive) as three separate 
characters, and partly a result of omitting taxa that 
linked groups differently in the Donoghue & Doyle 
analysis (e.g., Trimeniaceae, which tend to link 
Chloranthaceae with other Laurales: cf. Endress, 
1987). 

The apparent conflicts between the morpholog- 
ical and molecular data, especially regarding “root- 
ing” of the angiosperms, were discussed by Don- 
oghue & Doyle (1989b), who suggested that they 
may illustrate a general conclusion drawn by Hillis 
(1987): that conflicts between morphological and 
molecular results rarely reflect serious contradic- 
tions, but rather different levels of resolution of the 
two sorts of data in different parts of the tree. This 
seemed to be supported by the fact that exper 
ments with alternative trees in the seed plant stud- 
ies (Doyle & Donoghue, 1986, 1992) had shown 
that other positions of anthophytes within se 
plants, and of angiosperms within anthophytes, жыр 
only slightly less parsimonious than the most e 
simonious arrangements. For example, it was 0 y 
one step less parsimonious to link angiosperms 
Caytonia and/or glossopterids, thus brotas i 
the anthophytes, or to link anthophytes with Men 
eropsids, with Gnetales basal in anthophytee- f 
we called a neo-englerian arrangement. d 
in our morphological analysis of angiosperms ( еді 
oghue & Doyle, 1989a), we had found trees es 
angiosperms rooted among paleoherbs that y 
only one step less parsimonious than those xa 
near Magnoliales. This may therefore be A = 
where the morphological data are aE its 
they favor a woody magnoliid prototype» gene 
slightly—whereas the molecular data point $ E 
ly in one direction, toward à paleoherb rool! ; 
and therefore provide better evidence on 


data 
ships. In other cases, it may be the molecular 


with 
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that are ambiguous and the morphological data 
that are clearcut. Thus the two sorts of data may 
be more complementary than contradictory. 

In the present study, we have attempted to 
address these issues by comparing results derived 
from morphological and molecular data for the 
same set of seed plant taxa, probing the strengths 
and weaknesses of the results with methods such 
as bootstrap and decay analysis, and analyzing a 
combined data set. The question of analyzing mor- 
phological and molecular data separately and com- 
paring the results or combining them at the outset 
isa topic of ongoing debate (Kluge, 1989; Barrett 
et al., 1991; Donoghue & Sanderson, 1992; Bull 
et al., 1993; de Queiroz, 1993). One argument 
against combining data sets is that the greater 
number of molecular characters will simply over- 
whelm the morphological characters. However, this 
does not necessarily hold: if the molecular results 
are poorly resolved, as they often are, even a small 
number of morphological characters can have a 
decisive effect (Donoghue & Sanderson, 1992). In 
any case, it is possible both to analyze data sets 
separately and to combine them, and our results 
27 that this approach can give instructive ге- 


This study is not intended to be a comprehensive 
cao of "morphological" versus **molecu- 
е; data on this topic. Y. Suh (pers. comm.) has 
Uem different results from another part of the 

subunit of rDNA, which roots the angiosperms 

n a clade including most Magnoliales and 
е and other angiosperms. Equally different 
iid Ie the aquatic genus Ceratophyllum basal 
asec,’ angiosperms divided into tricolpate- 
ecu eudicots and monosulcate magnoliids and 
(les py have been obtained from rbcL sequences 
1993) PA 1991; Chase et al., 1993; Qiu et al., 

e : nalyses of shorter rRNA sequences by 

usky et al. (1991) and a smaller rbcL data set 
riii al. (1992) have given different trees, 
£p. I y with gymnosperms as a monophyletic 
e "i nstead of considering all these data sets, 

ve chosen to address the two that we know 
“st, although we will mention briefly some prelim- 
mary analyses including rbcL. 


ComBininc Pn 


EVIOUS MORPHOLOGICAL ANALYSES 
OF SEED Pran 


TS AND ANGIOSPERMS 


is eal concern that we wish to address first 
ing рна that previous inferences concern- 
by circu "un of angiosperms were compromised 
М t reasoning. Perhaps the position of an- 

rms in the seed plant analysis (Doyle & Don- 


oghue, 1986) was a function of the assumption 
that the first angiosperms were like Magnoliales 
and Winteraceae. Perhaps then the basal position 
of Magnoliales in the angiosperm analysis (Dono- 
ghue & Doyle, 1989a) was a consequence of an 
incorrect identification of outgroups based on this 
initial assumption. For the sake of rapid progress, 
we split the seed plant and angiosperm problems 
in two. This seemed reasonable at the time: there 
was already strong evidence that the angiosperms 
were monophyletic, and there was some consensus 
on basic states within the group. However, we 
realized from the beginning that it would eventually 
be necessary to carry out additional analyses de- 
signed to resolve simultaneously relationships at 
the point where the two analyses intersect. For 
example, we noted that Chloranthaceae share many 
features with Gnetales, and that angiosperms might 
therefore be directly associated with Gnetales if 
Chloranthaceae were assumed to be primitive. 

It should be recognized that conclusions derived 

from this sort of two-step procedure are not nec- 
essarily incorrect: there might be only one most 
parsimonious position for angiosperms no matter 
what internal relationships or basal states in an- 
giosperms are assumed. It should also be noted that 
the problems are not unique to our study: ours is 
simply one example of a general method that Mish- 
ler (1994) calls compartmentalization. This method 
was also used within our angiosperm analysis (Don- 
oghue & Doyle, 19893), in which we did a prelim- 
inary analysis of Laurales to determine basic states 
in a derived subgroup referred to as “соге Laura- 
les." 
These problems were recognized and addressed 
in a series of experiments reported by Doyle & 
Donoghue (1990), which combined nine angio- 
sperm taxa with the 17 nonangiospermous groups 
used in Doyle & Donoghue (1992). Here we pre- 
sent an updated version of these experiments, which 
illustrate the nature of the problem and the poten- 
tial value of considering both morphological and 
molecular evidence. 


TAXA, CHARACTERS, AND ANALYSES 


The revised seed plant and angiosperm matrix, 
henceforth designated the nine-angiosperm anal- 
ysis, is presented in the Appendix (Table 1). The 
angiosperm taxa were selected to represent the 
major clades found by Donoghue & Doyle (19892) 
in trees rooted both near Magnoliales and among 
paleoherbs, which necessitated dividing the pa- 
leoherbs into four groups. In five cases these clades 
are represented by individual taxa used in Dono- 
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ghue & Doyle (1989a), in four by composite taxa: 
Magnoliales, based on Degeneria, Myristicaceae, 
Annonaceae, and Magnoliaceae (“соге Magnoli- 
ales” of Donoghue & Doyle, 1989a); Piperales 
(Piperaceae, Saururaceae); Nymphaeales (Nym- 
phaeaceae, Cabombaceae); and eudicots (Ranun- 
culidae, Nelumbo, Trochodendrales, Hamameli- 
dales). These are scored in terms of estimated 
ancestral states for the whole taxon (theoretical 
and practical problems in this procedure are dis- 
cussed further below). Thus, where the taxon con- 
sists of two taxa included in the previous analysis 
and these have different stages, the group was 
scored as uncertain (e.g., 0/1). In the case of 
binary characters, where uncertain (0/1) and un- 
known (?) are equivalent in tree construction, we 
usually distinguished between the two scorings in 
the matrix to indicate the nature of the uncertainty, 
but we did not try to weed out all cases where “?” 
had been used for uncertainty in the previous anal- 
yses. In Magnoliales, basic states were estimated 
based on the previous result that Degeneria is the 
sister group of the other three taxa, which them- 
selves form an unresolved trichotomy. Eudicots 
were assumed to consist of two sister clades, Ranun- 
culidae plus Nelumbo and Trochodendrales plus 
Hamamelidales. Laurales were represented by Aus- 
trobaileya, which was at or near the base of the 
order in our previous trees, and Chloranthaceae, 
which are of special interest because they have 
been widely discussed as possible primitive angio- 
sperms. Use of Austrobaileya to represent Laur- 
ales might be questioned, since it lacks many fea- 
tures commonly associated with the order. However, 
our results bear out its use in this way, since Aus- 
trobaileya is associated with Chloranthaceae in the 
trees obtained, as it was with the larger data set. 
Trees were rooted by including a taxon based on 
Devonian '"*progymnosperms" (Aneurophyton, Ar- 
chaeopteris) as outgroup. 

Characters are primarily a combination of those 
used in the seed plant analysis of Doyle & Don- 
oghue (1992) and those in the angiosperm analysis 
of Donoghue & Doyle (1989a) that are potentially 
informative for the taxa under consideration. To 
these we added several apomorphies that poten- 
tially hold angiosperms together as a monophyletic 
group, either as new characters or as additional 
states of existing seed plant characters (e.g., three- 
nucleate microgametophyte, complete loss of the 
megaspore wall). In some cases, features that vary 
within angiosperms are expressed as additional states 
of characters recognized in seed plants as a whole 
(e.g., palmate leaf venation). For angiosperm char- 
acters in which there are no clearly comparable 


states in other seed plants, we scored the latter as 
unknown; this is especially true of floral characters, 
where we hoped to avoid biasing the results by 
assuming questionable homologies of parts between 
angiosperms and other groups. These characters 
also form the basis for the morphological analysis 
of Recent taxa presented below as a counterpart 
of the rRNA analysis; for simplicity we retained 
nine states that were potentially informative in that 
matrix but autapomorphic in the present one (e.g., 
tetracytic stomates in Piperales, tricolpate pollen 
in eudicots). 

Multistate characters were unordered, except 
for two easily ordered quantitative characters (pol- 
len size, megaspore wall thickness). Two multistate 
characters involving exine structure deserve spe- 
cial consideration, since our decision not to order 
them had a significant impact on the results pre- 
sented below. 

In the angiosperm study (Donoghue & Doyle, 
1989a), we recognized an infratectal structure 
character with granular and columellar states and 
assumed that granular was primitive, based оп out- 
group comparison with other anthophytes. In the 
revised seed plant study (Doyle & Donoghue, 1992), 
we recognized spongy-alveolar, honeycomb-alveo- 
lar, and granular states and scored angiosperms аз 
granular, assuming that columellar structure 
evolved within the group. In combining the two 
data sets, our previous polarization of states within 
angiosperms could be preserved by ordering the 
character (spongy-honeycomb-granular-columel- 
lar). By placing two steps between alveolar and 
columellar, this ordering would bias against trees 
in which angiosperms are linked with alveolar out- 
groups (e.g., Caytonia, glossopterids), eee 
groups are basal in angiosperms, and granular 
structure in groups like Magnoliales is à peed 
gence with Bennettitales, Pentoxylon, and Gne 
tales. This scenario involves four steps if the exin? 
character is ordered but three if it is not- ath 
bias might be defended based on the coexisten 4 
of columellar and granular structure (and end 
tional states) within angiosperms (Walker, l 
Le Thomas, 1980-1981) and of alveolar and qe 
ular structure within conifers (Van Campo & күс 
gardon, 1973). However, it seems unwarrant ха 
assume that a direct transition from are 
columellar could not have occurred in other es 
Trees with angiosperms linked with Cayton/ 

in terms 
only one step longer than the shortest trees xn 
of the Doyle & Donoghue (1986, 1992) data *5/ 
and a transition from alveolar to granular 18 те 
tainly conceivable on structural grounds ыы 
duction of the side-walls of the alveolae, 
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the junctions between the walls as columellae: E. 
Masure, pers. comm.). 

Endexine structure was not included as a char- 
acter in the seed plant analyses (Doyle & Dono- 
ghue, 1986, 1992), because all groups except an- 
giosperms have a uniformly laminated endexine, 
making the character uninformative. Of the two 
states in angiosperms, endexine present or absent 
in the extra-apertural areas, absence of endexine 
was assumed to be ancestral (Donoghue & Doyle, 
1989a). This was based on the hypothesis that the 
laminated endexine of other seed plants is homol- 
ogous with the footlayer of angiosperms, since both 
develop from similar tangential lamellae, and that 
the endexine of angiosperms, which is nonlami- 
nated except sometimes under the apertures, is a 
new layer (Zavada, 1984). Again, the original po- 
larization could be preserved in the combined data 
set by ordering the character (endexine laminated- 
absent-nonlaminated). However, because the ho- 
mologies involved are rather speculative (cf. Ga- 
barayeva, 1991), it seems preferable to treat the 
three states as unordered. 

We also made a few substantive changes based 
on new data, such as recognition that pollen of 
Piperaceae (Piperales) has supratectal spinules and 
а sculptured sulcus (Bornstein, 1989), Cabomba- 
ceae (Nymphaeales) have a columellar exine struc- 
ture (Osborn et al., 1991), Myristicaceae (Mag- 
noliales) have both S and PI type sieve-tube plastids, 
ànd Aristolochiaceae have basically PII type plas- 
lids (as seen in Saruma and Asarum), like mono- 
tots (Behnke, 1988). We added one new potential 
‘ynapomorphy of angiosperms and Gnetales, dou- 
i fertilization, in the sense of regular fusion of 

* second sperm nucleus with a second megaga- 
metophyte nucleus. This has been confirmed in 
s Se and seems independent of (although prob- 
feat, oe pd for) the uniquely angiospermous 
1999. vom osperm formation (Friedman, 1990, 

+ Donoghue & Scheiner, 1992). 

the i Sets of analyses were performed using 
Gal ы nine-angiosperm matrix. In one set, we 
Yzed the matrix with each of the nine angio- 
‘Perm groups substituted individually for angio- 
ene removing the eight remaining taxa, we 
Ше characters that then became autapomor- 
Sian жы often desirable to remove such сһаг- 
oes ause they are uninformative and distort 
"pu ез of homoplasy, but this is not a problem 
with > case, where we are not concerned 
fiir a levels of homoplasy. However, we did 
hA variant characters with MacClade (Mad- 
tmi Maddison, 1992). All these data sets were 
yzed with PAUP (version 3.0L, Swofford, 


1991). Since there are too many taxa for branch- 
and-bound analysis, which guarantees finding all 
most parsimonious trees, we used the heuristic 
search algorithm, with 10 replicates of stepwise 
random addition of taxa and TBR branch swapping. 
This increases the probability of finding most par- 
simonious trees that belong to different “ізіапав” 
(Maddison, 1991), which were in fact found in 
some experiments (see results). Alternative ar- 
rangements were investigated using MacClade and 
the constraints option in PAUP. 


SINGLE-ANGIOSPERM ANALYSES 


Experiments with single angiosperm taxa re- 
sulted in several different positions of the angio- 
sperms relative to other anthophytes (Pentoxylon, 
Bennettitales, Gnetales), and of anthophytes within 
seed plants. 

As expected, when Magnoliales are substituted 
for angiosperms (Fig. 1), angiosperms are the sister 
group of other anthophytes, and anthophytes are 
associated with one or another combination of Cay- 
tonia, glossopterids, and corystosperms, as in Doyle 
& Donoghue (1986, 1992). Relationships among 
other groups also parallel those found by Doyle & 
Donoghue (1992), including some with the ar- 
rangement of extant cycads, Ginkgo, and conifers 
found by Loconte & Stevenson (1990). As shown 
in Figure 1, one of the characters that supports 
this position of Magnoliales is the presence of gran- 
ular exine structure. 

In contrast, when Winteraceae, Austrobaileya, 
eudicots, Aristolochiaceae, or monocots are sub- 
stituted for angiosperms, the angiosperm taxon 
connects with Caytonia or Caytonia plus glossop- 
terids. This is also true of most trees found when 
Nymphaeales are substituted for angiosperms (Fig. 
2), although in two such trees the Bennettitales- 
Pentoxylon-Gnetales clade, Nymphaeales, and 
Caytonia form a paraphyletic group at the base 
of the platysperms (a stratigraphically very un- 
parsimonious arrangement). These trees break up 
the anthophytes, although the Bennettitales~Pent- 
oxylon-Gnetales clade is still the next-closest group 
to angiosperms. This result might be questioned 
because it implicitly assumes that Caytonia and/ 
or glossopterids had angiosperm and gnetalian states 
for several characters that were scored as unknown 
because they are not preserved or not yet estab- 
lished in fossils, such as lignin chemistry, a tunica 
layer in the apical meristem, and siphonogamy (the 
same is also true for Bennettitales and Pentoxylon 
in trees of the sort shown in Fig. 1). One reason 
for the new position of these angiosperm taxa is 
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Treelength: 132 


65 Chars. E honeycomb 
ШЕШ granular 
FiGURE 1. Representative most parsimonious tree (of 


94) found when “соге” Magnoliales (MAGN) аге substi- 
tuted for angiosperms as a whole in the nine-angiosperm 
analysis of extant and fossil seed plants. Shading of branch- 
es shows distribution of the exine structure character, 
which tends to link Magnoliales with other anthophytes 
(Pentoxylon, Bennettitales, Gnetales). PROG = *'progym- 
nosperms"; ELKI = Elkinsia; MEDU = Medullosaceae; 
CALL = Callistophyton; CONI = Coniferales; GINK = 
Ginkgoales; CORD = Cordaitales; PELT = Peltasper- 
mum; CYCA = Cycadales; CORY = Corystospermaceae; 
GLOS = Glossopteridales; CAYT = Caytonia; PENT = 
Pentoxylon; BENN = Bennettitales; ЕРНЕ = Ephedra; 
WELW = Welwitschia; GNET = Gnetum. 


presumably that they have columellar rather than 
granular exine structure (or both columellar and 
granular structure in the case of Nymphaeales, 
which were therefore scored as uncertain), so they 
аге not as strongly "attracted" to other granular 
anthophytes as Magnoliales are. The exine struc- 
ture character would have tended to associate an- 
giosperms with other anthophytes if it had been 
ordered, with columellar implicitly derived from 
granular, but it does not when it is treated as 
unordered. Competing characters attracting angio- 
sperms in general to Caytonia are reticulate ve- 
nation, flat guard cells, and anatropous cupules 
(scored like anatropous bitegmic ovules). 

As anticipated, in most of the trees found when 
Chloranthaceae are substituted for angiosperms (Fig. 
3a), Chloranthaceae are linked directly with Gne- 
tales, with which they share such features as op- 
posite leaves, two-trace nodes, spicate inflores- 
cences (scored as compound strobili), and 
orthotropous ovules. In these trees, the position of 
anthophytes is highly unstable: they may be the 
sister group of other platysperms or variously as- 
sociated with cycads, glossopterids, Caytonia, or 
a corystosperm-glossopterid- Ca ytonia clade. The 
exceptions are neo-englerian trees in which Chlo- 
ranthaceae are linked with Pentoxylon and Ben- 





Treelength: 130 
66 Chars. 
= ЕЕЗ equivocal 


Ғісиве 2. Representative most parsimonious tree (of 
31) found when Nymphaeales (NYMP) are substituted 
for angiosperms, showing distribution of the exine struc: 
ture character. Generally similar trees (some with angio- 
sperms linked with Caytonia alone) are also found when 
Winteraceae, Austrobaileya, eudicots, Aristolochiaceae, 
and monocots are substituted for angiosperms. Arrows 
indicate possible exine states on branches where the state 
is equivocal, Other abbreviations as in Figure it 


nettitales (Fig. 3b), which maintain the putative 
homologies between Chloranthaceae and Gnetales, 
such as compound strobili, as symplesiomorphies. 
Anthophytes are not broken up in any of these 
trees. ; 

The most varied trees are found when Piperales 
are substituted for angiosperms. As when Chloran- 
thaceae are the single angiosperm group, these 
include neo-englerian trees where Piperales are 
linked with Pentoxylon and Bennettitales and cé: 
where anthophytes are the sister group of ot е 
platysperms. However, іп trees of the latter 1 
Piperales are not linked with Gnetales but un 
with Pentoxylon and Bennettitales (Fig. 4). я 
is presumably because Piperales have eunte 
orthotropous and bitegmic ovules that allow ыл? 
to be nested between Bennettitales and Gnet - 
but not features like opposite leaves and derit 
nodes that unite Gnetales and Chloranthecist 
However, Piperales are linked with Gnetales a 
few stratigraphically unparsimonious pu e 
anthophytes form a paraphyletic group а! 
of platysperms. Finally, there are ae 
in which Piperales are basal in anthophyt Wer 
anthophytes are nested among Caytonta, ae of 
terids, and corystosperms, analogous 10 when 
Doyle & Donoghue (1992) and trees fo 
Magnoliales are substituted for angios 
1), and a few trees where Piperales are 


(Fig 
pw 
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FICURE3. Two representative most parsimonious trees 


Lori found when Chloranthaceae (CHLO) are substi- 
мах (a angiosperms, showing distribution of the phyl. 
ғ: у а)апа strobilus (b) characters; (b) is а neo-englerian 

» with anthophytes nested among “‘coniferopsids” 


он. conifers, cordaites). Other abbreviations as in 


nima ES Caytonia, thus breaking up the an- 
pP resumably, Piperales have fewer char- 
Ras support any one of these arrangements 
another, 

in pend of angiosperm relationships seen 
ib periments is not evident when extant 
Castel, 2 = ALI too since in almost all trees 
Sorina (i е the closest living relatives of angio- 
Қы а exceptions аге the few trees in which 
group), д and Gnetales form a paraphyletic 
rationali quus are a function of different 
"we - tween angiosperms and various fos- 
very differe ы; the different trees would have 
ater =a implications for basic states and char- 
astra: ion within the angiosperms, and this 
extant 2 s limitations of trees based on 
scenarios ie one in evaluation of evolutionary 
groups е cause of the gaps between extant 

: For example, trees in which angiosperms 
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FIGURE 4. Representative most parsimonious tree (of 
155) found when Piperales (PIPE) are substituted for 
angiosperms, showing distribution of the cupule character. 
Other equally parsimonious trees are more comparable to 
those in Figures 1, 2, and 3b. Other abbreviations as in 


Figure 1. 


are linked with Caytonia and/or glossopterids im- 
ply that flowerlike reproductive structures origi- 
nated independently in angiosperms and the clade 
consisting of Bennettitales, Pentoxylon, and Gne- 
tales. Trees in which angiosperms are linked with 
Gnetales, and/or neo-englerian trees in which an- 
thophytes are linked with coniferopsids and Gne- 
tales are basal, imply that angiosperm flowers and 
floral parts, especially carpels containing several 
ovules with two integuments, were elaborated from 
simpler structures like those of Gnetales, or derived 
by aggregation of several gnetalian “flowers” (i.e., 
a pseudanthial interpretation). 

These results bear out the concern that the 
position of angiosperms in previous analyses (Doyle 
& Donoghue, 1986, 1992) may have been incor- 
rect because of initial assumptions about basal states 
in angiosperms. They also raise the possibility that 
angiosperms are polyphyletic, with different “ап- 
giosperm" groups related to different "gymno- 
sperms.” However, neither conclusion necessarily 
follows, since they depend on whether and how the 
various angiosperm groups link up with each other. 


NINE-ANGIOSPERM ANALYSES 


To assess the possibilities just raised, we included 
all nine angiosperm groups and analyzed the re- 
sulting matrix (35 replicates, stepwise random ad- 
dition of taxa, TBR branch swapping). This analysis 
yielded 11 most parsimonious trees of 192 steps, 
which differ only in arrangements within a clade 
consisting of Callistophyton, coniferopsids, cory- 
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FiGURE 5. Representative most parsimonious tree (of 11) found when all nine angiosperm taxa are saa р 
the analysis, showing distribution of the exine structure character. PIPE = Piperales; ARIS = At i» 
= Nymphaeales; MONO = monocots; MAGN = “core” Magnoliales; WINT = Winteraceae; EUDI т ещ u s 
groups with tricolpate and derived pollen); AUST = Austrobaileya; CHLO = Chloranthaceae; other abbrevia 


in Figure 1. 


stosperms, Peltaspermum, and cycads (e.g., Fig. 
5). In these trees, angiosperms form a monophy- 
letic group, and Gnetales are their closest living 
relatives. However, in contrast to trees of Doyle 
& Donoghue (1986), where angiosperms are the 
sister group of other anthophytes, angiosperms are 
linked with Caytonia and glossopterids, as when 
Winteraceae, Austrobaileya, eudicots, and some 
paleoherbs were substituted for angiosperms (Fig. 
2). Furthermore, the arrangement within angio- 
sperms differs from any seen previously: they split 
into one clade consisting of paleoherbs (including 
monocots) and another consisting of woody mag- 
noliids and eudicots. As with trees where single 
angiosperm taxa were associated with Caytonia 
and/or glossopterids, this rearrangement is pre- 
sumably influenced by treatment of the exine struc- 
ture and endexine characters as unordered. This 
change in character analysis weakens the tendency 
of angiosperms to associate with the granular Ben- 
nettitales- Pentoxylon-Gnetales clade, with gran- 
ular Magnoliales attached between the latter and 
columellar angiosperms. 

Although these results cast doubt on previous 
inferences regarding the origin of angiosperms, the 
change from the previous situation is less radical 
than it seems. Essentially the same alternative re- 
lationships recognized as almost equally parsimo- 


nious by Doyle & Donoghue (1986, 1992) and 
Donoghue & Doyle (1989a) are seen among еу 
step less parsimonious (“опе-оЖҒ”) trees (193 кере 
(Fig. 6). Some of these trees show the sorts 
relationships that were most parsimonious in ! 
previous analyses, in which anthophytes are ае 
sociated with Caytonia and glossopterids, angi” 
sperms are the sister group of other eae oh 
and Magnoliales are basal in angiosperms (Fig. 4 
As іп the previous analyses, other one-off ке : 
of the neo-englerian type (Fig. 6b), in whic еі 
thophytes are linked with coniferopsids, and А 
tales, with linear leaves and simple sporophylls, де 
basal and relatively plesiomorphic іп wenn 5 
In other 193-step trees, Nymphaeales аге е 
angiosperms, as inferred from rRNA data (Ha e 
& Zimmer, 1992); in fact, some of these have ds 
same arrangement of paleoherb taxa irn uii 
rRNA analyses presented below, with Piper pe 
above Nymphaeales and monocots linked us реуі 
tolochiaceae (Fig. 6c). Given the minim игу 
ences in parsimony, it would be gram м 
conclude that the tree in Figure 5 een 
strongly preferred over those in Did 4 
(19894) or Doyle & Donoghue (1992). р 
the shift away from a basal position of — 
may represent real progress, reflecting бөгесін 
the removal of previous biases іп ехіле cha 
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A es 6. Representative one step less parsimonious (one-off) trees found with the nine-angiosperm data set: (a) 
angiosperms the sister group of other anthophytes and Magnoliales basal in angiosperms, as in Doyle & Donoghue 


(1 ) 
ed 1992) and Donoghue & Doyle (1989), showing distribution of the exine structure character; (b) neo-englerian 
ang with anthophytes nested among "coniferopsids," showing distribution of the sporophyll character; (c) with 

Naberms linked with Caytonia and glossopterids and Nymphaeales basal in angiosperms, as in trees based on 


"RN; . . М d 
A data, showing distribution of the exine structure character. Abbreviations as in Figures 1 and 5. 
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FicURE 7. Representative tree (of 1675) found when angiosperms and Gnetales are forced together and aci 
ranthaceae are forced to the base of angiosperms, showing distribution of the phyllotaxy character. This tree 
steps less parsimonious than the shortest trees. Abbreviations as in Figures 1 and 5. 


it may be as much a function of the smaller number 
of taxa and the omission of connecting groups. 

Despite these ambiguities, experiments with al- 
ternative topologies do provide some stronger re- 
sults. The most interesting experiments concern 
proposed links between Gnetales and Chlorantha- 
ceae. Trees with angiosperms forced together with 
Gnetales using the constraints option in PAUP are 
only two steps longer than the shortest trees, but 
Chloranthaceae are not basal in angiosperms, as 
might be expected from their gnetalian features; 
instead, angiosperms are arranged as in the most 
parsimonious trees (Fig. 5). The shortest trees with 
Chloranthaceae basal in angiosperms and angio- 
sperms linked with Gnetales (Fig. 7) are six steps 
longer than the shortest trees. Similarly, in the neo- 
englerian trees (Fig. 6b), anthophytes are not ar- 
ranged in such a way that the similarities between 
Gnetales and Chloranthaceae are homologous; in- 
stead, angiosperms are linked with Bennettitales, 
and Magnoliales are basal in angiosperms. 

The reason that trees of the type in Figure 7 
are so unparsimonious, even though Chlorantha- 
ceae are associated with Gnetales when they are 
the only angiosperms in the analysis, is presumably 
that Chloranthaceae are “screened off” from hav- 
ing an influence on the position of angiosperms by 
being linked with Austrobaileya and other groups, 
based on such features as laterocytic stomates and 


globose pollen with reticulate sculpture and a ver- 
rucate sulcus. Similarly, Piperales are nested among 
other paleoherbs based on herbaceous habit, " 
tichous leaves, palmate venation, and globose pot- 
len. Of the potential homologies of mace 
and Gnetales, opposite leaves and two-trace П es 
are unequivocally primitive in angiosperms In es 
like Figure 7, since they also occur in eo A 
baileya. However, compound strobili would a а 
to be lost below Austrobaileya, so it ien 
parsimonious to assume that this feature € 
independently in Chloranthaceae and Gnet = 
Other “gnetalian”’ features of Cention E 
problematical morphologically and were ves 
scored as unknown (e.g., whorled тісгоѕрогор ^ 
one ovule per carpel). As a result, they are 
sistent with a link between Chiari Y: 
Gnetales when Chloranthaceae alone = ваг: 
tuted for angiosperms, but they do not "at re 
Chloranthaceae to the base of and al 
all nine angiosperm taxa are included. Fin уе 
though Chloranthaceae are like Gnetales we z 
respects, Magnoliales are more like rk өлі» 
others (e.g., leaves without chloranthoid tee gm 
shaped pollen; tectate, granular exine t oe 5 
The strongest conclusion of these en 
that both angiosperms and Gnetales MS mo they 
letic groups. In the tree shown in Figure > - 


are unambiguously supported by nine 8 
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characters, respectively. Four of the angiosperm 
synapomorphies are universal in the group (or rep- 
resented by clearly related states) and unknown 
elsewhere (two pairs of pollen sacs, endothecium, 
stigmatic pollen germination, loss of megaspore 
wall). Of the five others, vessels and several vein 
orders arose independently in Gnetales, scalariform 
secondary xylem in Bennettitales and Pentoxylon; 
columellae and nonlaminated endexine are not uni- 
versal in angiosperms but are basic with this ar- 
rangement. Four other features that are known 
only in angiosperms are equivocal as angiosperm 
synapomorphies because the corresponding char- 
acters are unknown in Caytonia and glossopterids 
(companion cells, three-nucleate microgameto- 
phyte, eight-nucleate megagametophyte, endo- 
sperm formation). To evaluate the strength of this 
result, we removed angiosperm apomorphies to see 
at what point the group would break up. Remark- 
ably, even when we removed all eight features that 
are known only in angiosperms, angiosperms still 
stayed together as a clade. Apparently there are 
enough overlapping similarities within angiosperms, 
such as trilacunar nodes, columellar exine struc- 
ture, ой cells іп most magnoliids, and palmate ve- 
nation in paleoherbs and eudicots, to hold them 
logether. 

To some, this result may seem trivial, but even 
iie authors have expressed the opinion 
dei зна monophyly is а pernicious dogma 
2s s held back progress in understanding the 
етенсе (e.g., Hughes & McDougall, 
Cann assiloy, 1991). Similarly, the view that 
Aes chm polyphyletic is still frequently en- 
idis TR | (e.g., Gifford & Foster, 1989). Our re- 
es € y that. the assumption that angiosperms 
фа "Piin oi is not an obstacle to progress in 

» but incorrect assumptions concerning the 
рад of the first angiosperms (an inappro- 
“search image") might well be. 
tat ese results underline a general problem іп 
ei Бои that are separated from 
discum As by “опр branches”: the more 
i its mon E unite the group, the more certain 
itt posit op yletic status, but the less certain is 
б е This is because spurious convergences, 
Seamed changes leading to uninterpretability 
relations; on the long branch may obscure true 
bis uo (Felsenstein, 1978). This effect has 
t often stressed in connection with mo- 
raters, where there are only three al- 
x ае to which а base at any position can 
ee t it may also apply to morphological 
nts (cf. Wake, 1991). For the same rea- 
» Tooting of the group may be ambiguous, be- 


cause the closest outgroups are so distant that they 
provide little signal" as to which subgroups are 
basal (Wheeler, 1990). This problem is reflected 
in the large number of angiosperm characters that 
could not be polarized by outgroup comparison in 
Donoghue & Doyle (1989a), or could not be scored 
outside of angiosperms in the present analysis, and 
it is magnified by the large number of missing 
characters in fossils. As discussed in greater detail 
elsewhere (Doyle & Donoghue, 1993), the problem 
might be solved by discovery of fossils on the long 
branch leading to angiosperms (i.e., non-angio- 
spermous angiophytes, or stem angiophytes, in the 
terminology of Doyle & Donoghue, 1993). How- 
ever, although there are a few candidates, like the 
Late Triassic Crinopolles pollen group described by 
Cornet (19892), which has angiospermlike retic- 
ulate sculpture and columellae but a gymnosperm- 
like endexine (cf. Doyle & Hotton, 1991), there 
are still no fossils with angiosperm states in some 
characters and more plesiomorphic states in others 
that can be placed with certainty on the angiosperm 
stem lineage. 


MORPHOLOGICAL AND rRNA ANALYSES OF 
ExTANT GROUPS 


The relationship of these results to those ob- 
tained from rRNA sequences is addressed more 
directly by our analyses of morphological and rRNA 
data from the same taxa. Since position and rooting 
of the angiosperms are ambiguous with the mor- 
phological data set just presented, and since still 
other results have been obtained with other inter- 
pretations of angiosperm characters (Loconte & 
Stevenson, 1991; Taylor & Hickey, 1992), one 
of the motivations for this study was to determine 
what if any additional insights into these questions 
can be extracted from rRNA data. Although mo- 
lecular data have the disadvantage of being avail- 
able only from living groups (except for a few recent 
fossils: e.g., Golenberg et al., 1990), whereas mor- 
phological data may exist for key fossil taxa that 
attach to the long branches separating extant groups 
(Donoghue & Doyle, 1989b; Donoghue et al., 
1989), molecular data have the advantage of being 
independent of the seemingly endless controversies 
over interpretation of the morphological homolo- 
gies of angiosperm structures. 


TAXA, CHARACTERS, AND ANALYSES 


The starting point for our rRNA analyses was 
a successor to the 60-taxon data set of Hamby & 
Zimmer (1992), enlarged to include 71 taxa. Com- 
plete sequences are available from GenBank (Ac- 
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cession Nos. M81965-M82800) and the NMNH 
Gopher Server under *LMS." Our morphological 
data set is derived from the nine-angiosperm anal- 
ysis described above, modified by inclusion of some- 
what different angiosperm taxa and removal of 
fossil groups (Appendix). 

We encountered a variety of problems in ob- 
taining comparable taxa for the morphological and 
rRNA analyses and combining the two data sets. 
This required many compromises and approxi- 
mations. All of these involve a certain risk of error, 
but we feel they are unavoidable if progress is to 
be made at this point (cf. Maddison & Maddison, 
1992). The important thing is to spell out the 
assumptions involved so that they can be scruti- 
nized and tested in future analyses. 

In reducing the morphological and molecular 
data sets to a comparable set of taxa, our goal was 
to include an adequate sampling of critical taxa (as 
judged from results of previous studies) while keep- 
ing the number of taxa small enough for the more 
time-consuming analyses. We had to omit taxa for 
which rRNA data are still lacking; this is often 
unfortunate, since current evidence suggests that 
some such groups constitute important links. Ex- 
amples are Lactoris, which recent authors have 
linked with Piperales (Carlquist, 1990) or Aristo- 
lochiaceae (Donoghue & Doyle, 19892; Qiu et al., 
1993), and Austrobaileya, Trimeniaceae, and Am- 
borella, which may help to tie together Laurales 
and strengthen the position of Chloranthaceae in 
this group. We did not include Ceratophyllum (also 
omitted by Donoghue & Doyle, 19892), because 
its position was highly unstable in the complete 
rRNA analyses (Hamby & Zimmer, 1992) and so 
many of its morphological characters are difficult 
to interpret. Although Ceratophyllum occupies a 
key position at the base of the angiosperms in rbcL 
analyses of seed plants as a whole (Les et al., 1991; 
Chase et al., 1993), its-position is unstable in un- 
rooted rbcL analyses of angiosperms alone (Qiu et 
al., 1993), suggesting that its basal position may 
be an artifact of long branch attraction (cf. Don- 
oghue, 1994). 

Similarly, we did not include any non-seed plants 
as outgroups, since outgroups in the original data 
sets were different. In the morphological analyses 
of extant seed plants (Doyle & Donoghue, 1987, 
1992), we assumed that ferns and Equisetum are 
closer to extant seed plants and lycopsids more 
distant, but the outgroups included in the rRNA 
analysis were Equisetum and Psilotum. In addition, 
all extant outgroups are very distant from seed 
plants; the really appropriate outgroups are De- 
vonian "progymnosperms" and Carboniferous "seed 


ferns." There is therefore reason to fear that any 
rooting obtained from extant outgroups might be 
an artifact of spurious long branch attraction 
(Wheeler, 1990; Maddison et al., 1992). The re- 
sulting trees are therefore unrooted networks, with 
the root arbitrarily placed along the branch leading 
to cycads. However, it will be seen that this pro- 
cedure does provide instructive contrasts, because 
trees derived from the morphological and rRNA 
data sets are topologically different. 

In combining original taxa into larger groups, 
we generally accepted clades that appeared in both 
the morphological analysis of Donoghue & Doyle 
(19892) and the consensus of most parsimonious 
and one-off rRNA trees, with a few exceptions 
motivated by a desire to test current hypotheses, 
Thus we retained Piperaceae and Saururaceae às 
separate taxa, even though they were strongly linked 
in the trees of Donoghue & Doyle (19892) and 
associated in some one-off rRNA trees, because 
they are separated in the tree of Taylor & Hickey 
(1992). Conversely, we combined Hedycarya 
(Monimiaceae) and Persea (Lauraceae) as “core 
Laurales” in the rRNA analysis, even though they 
are not associated in all most parsimonious rRNA 
trees, because their relationship is strongly Sup- 
ported and uncontroversial on morphological 

ounds. е 
д Despite our efforts, taxa іп the morphological 
and rRNA data sets are not perfectly comparable, 
except perhaps when they are monotypic (е. 
Ginkgo, Welwitschia). In general, clades ш 
morphological data set аге represented in the ек 
lecular data set by а few species that show 7 
part of the variation in the whole clade (е.в ан 
noliales by Asimina, Magnolia, and Lirioden is 
core Laurales by Persea and Hedycarya), or та 
single “exemplar” species (Chloranthaces®, 7 
Chloranthus; Ranunculidae by Камаа 
Trochodendrales by Trochodendron). One so A 
would be to rescore taxa in the morphological 
set to correspond exactly to the species п kd 
analysis. However, we opted instead to score pe 
as in the earlier morphological sot ci 
assume that they are adequately represent 
rRNA analysis by the exemplars, since Lye 
unlikely that rescoring them would lead ы ; 
icantly different results іп the new ee Trace 
analysis. For example, Ranunculus 
dendron have most of the characters peni ac 
link Ranunculidae, Trochodendrales, and o 
dicots in the previous analysis cages 
sculptured aperture membranes, chloranthow” 
lack of oil cells), plus one that would Бе: в 
if Ranunculidae as a whole were co 
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gud о within that group (stamens with 
la ntiated filaments). 

in EU taxa that vary for characters included 
мо 1х, our goal was to obtain the best 
Xo ння states for the whole taxon. 
бажа thes ata set, we accepted internal rela- 
taxon rRNA a Pw consistent in the original 71- 
by. ourselves na Sigs: and in morphological analyses 
considerin xd others (Fig. 8a). For example, 
sumed mae = included in the analysis, we as- 
Brepols. is the sister group of Zamia and 
Арел pa Pinus is the sister group of 
are found in th -ryptomeria, since these relations 
morphological e whole rRNA analysis and in the 
(1988), са ts of Hart (1987), Crane 
соне lera (1990). In contrast, in 
claya is the z e rRNA data indicate that Bar- 
but Hr MER group of the remaining taxa, 
vie nes ‘Gata (Ito, 1987) imply that Ca- 
td the three сиру this position; therefore we treat- 
groups as a trichotomy. In estimating 


basic states for monocots, we accepted relationships 
among grasses that are consistent with rRNA and 
morphological analyses (Kellogg & Campbell, 
1987), treated Hosta, Sabal, and grasses as a 
trichotomy, and treated alismids, aroids, and the 
Hosta-Sabal-grass clade as a trichotomy. In the 
morphological data set, we inferred states for core 
Magnoliales based on the assumption that Dege- 
neria is the sister group of Myristicaceae, Anno- 
naceae, and Magnoliaceae (cf. above); in core 
Laurales, we assumed relationships shown in Figure 
8b, derived from an unpublished analysis of Laura- 
les used for the same purpose by Donoghue & 
Doyle (1989a). 

The use of parsimony in optimization of ances- 
tral states on trees is discussed by Swofford & 
Maddison (1987). As in the nine-angiosperm anal- 
ysis, when clades consisted of two taxa that varied 
at a given site, we coded them as uncertain. When 
there were three taxa, and the "outer" taxon and 
one of the two "inner" taxa had the same state, 
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we interpreted this state as ancestral; but when the 
outer taxon differed from both inner taxa, we scored 
the clade as uncertain. In more complicated cases, 
we used MacClade (Maddison & Maddison, 1992) 
to find the most parsimonious ancestral state, tak- 
ing into account all possible resolutions of trichot- 
omies and polychotomies. 

The practice of including variable taxa and scor- 
ing ancestral states as uncertain has been criticized 
by Nixon & Davis (1991). In the absence of ho- 
moplasy, the main effect of scoring taxa as un- 
certain is lowered resolution, but Nixon & Davis 
presented theoretical cases where uncertainties 
combined with homoplasy lead to incorrect trees, 
a danger discussed in greater depth by Maddison 
& Maddison (1992: 47-49). The alternative pro- 
posed by Nixon & Davis (1991) is to split up 
variable taxa into units that are monomorphic in 
terms of the characters used (a procedure with 
risks of its own: Maddison & Maddison, 1992; 
Donoghue, 1994). However, in practice Nixon et 
al. (1994) used smaller taxa as exemplars for larger 
clades. As illustrated graphically by the experi- 
ments described above, where we obtained widely 
varying trees when we substituted different sub- 
groups for angiosperms as a whole, this approach 
is ridden with as many implicit assumptions and 
potentials for error (e.g., that convergences and 
reversals in the exemplar will not affect its position) 
as attempts to reconstruct basal states (Donoghue, 
1994). Splitting up variable taxa into all potentially 
relevant monophyletic taxa might be the ideal so- 
lution, but if carried to its logical conclusion this 
would quickly lead to computational paralysis. In 
the meantime, we prefer to make explicit assump- 
tions about basal states of the sort described, while 
emphasizing that these assumptions can and should 
be tested in the future. 

Our morphological data set is presented in Table 
2 (Appendix). Characters are the same as those in 
the nine-angiosperm analysis; changes in the num- 
ber and/or definition of states as a result of removal 
of fossils are indicated in the character definitions. 
Deletion of characters that became uninformative 
after removal of fossil taxa left a total of 69 char- 
acters. This data set therefore parallels the extant 
seed plant matrix of Doyle & Donoghue (1992). 

The rRNA data set is presented in Table 3 
(Appendix). Characters from the 18S subunit are 
keyed to the corresponding positions in soy, char- 
acters from the 26S subunit to positions in rice. 
As a result of reducing the number of taxa from 
71 to 18 and eliminating characters that became 
uninformative, the number of characters was re- 
duced from 411 to 174, of which 167 are base 


substitutions and seven are insertion-deletion events 
(indels). We included the indels in the analysis 
because they are all of short length (five involve 
single nucleotides, one a dinucleotide, one a tet- 
ranucleotide: see Appendix) and because the flank- 
ing sequences leave no ambiguity as to their align- 
ment. 

All three data sets were analyzed with the heu- 
ristic algorithm in PAUP (Swofford, 1991), with 
100 replicates of stepwise random addition of taxa 
and TBR branch swapping. Alternative topologies 
were investigated with the constraints option in 
PAUP and with MacClade (Maddison & Maddison, 
1992). 

One method used to evaluate the relative strength 
of various results is bootstrap analysis (Felsenstein, 
1985). Characters are sampled randomly from the 
original data set with replacement and trees are 
calculated for the new data set, and this procedure 
is replicated many times. The original concept was 
that the frequency of bootstrap replicates in which 
a clade occurs is an estimate of its statistical sig- 
nificance: e.g., if two taxa are united in 95 out of 
100 replicates, their relationship can be accepted 
at a 95% confidence level. Whether bootstrap fre- 
quencies should be interpreted in this way 15 4 
matter of debate (Carpenter, 1992; Hillis & Bull, 
1993; Felsenstein & Kishino, 1993), but ће спі. 
icisms made do not call into question the value of 
bootstrap analysis as a means of evaluating the 
relative robustness of clades. Actually, simulation 
experiments suggest that the bootstrap errs 0n the 
side of being conservative; under many circum- 
stances, clades seen at bootstrap levels lower than 
the conventional limit of 95% are more eee 
than the bootstrap numbers would imply (Hillis 
Bull, 1993). In addition, bootstrap analysis ip 
be useful in uncovering possible alternative 16 А 

$ Р ^ T f “frequency 9 
tionships, as seen in the dot-plots o 

» i PAUP. Ifa 
occurrence” of groups provided by PA 
link seen іп the most parsimonious trees 15 317 d 
due to convergence, there should also be minority 
characters that reflect the true relationships t 
these should be sampled and amplified in 4 For 
tively high frequency of bootstrap ropie 
each data set, we did 1000 bootstrap MeS 
To increase the probability of finding pA 
monious trees in each replicate, we did 10h 
analyses with stepwise random addition of taxa 
TBR branch swapping. 

The other ied Ae employed to n 
bustness is decay analysis (Bremer, 1988; o 
ghue et al., 1992). This method da 
much longer trees have to be—how muc found 
mony has to be relaxed— before trees аге 
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Fi m c. z 
Ж сан Two of the four most parsimonious trees found in analysis of the extant morphological data set, showing 
г of unambiguous changes supporting each clade. Abbreviations as in Figure 8. 


ei a given clade breaks up. With PAUP, 
c x Vu ri by retaining all trees equal to 
pie than a given length, constructing a strict 
Mic of the resulting trees, and observing 
P clades remain in the consensus. Ideally, this 
one with a branch-and-bound algorithm, which 
cs e finding all trees of a given length; be- 
rel was not possible with 18 taxa, we used 
id tic search with 100 replicates of stepwise 
om addition and TBR branch swapping. 


R 
ESULTS OF MORPHOLOGICAL ANALYSES 


ete of the morphological data set yields 
ы parsimonious trees of 150 steps, two 
whether Mis Figure 9. The variations concern 
ME Ua are the sister group of the 
Vibe ^is lan taxa or nested among them. Figure 
MR. the number of characters that un- 
piene Y support clades; angiosperms are united 
Gite $ fifteen characters, Gnetales by eight, 
indication Ns groups by seven. This gives some 
mised; о Sn level of support, but it is potentially 
tribution d cause it says nothing about the dis- 
are uni У homoplasy—whether these changes 
This ae or duplicated elsewhere on the tree. 
бағы ы is addressed by the bootstrap and 
оу nalyses. The consistency index is 0.58, 
average for this number of taxa (Sanderson 
DUM, 1989); the retention index is 0.73. 
би, а show the arrangement of non-an- 
"s groups found in analyses of extant 


taxa alone by Loconte & Stevenson (1990) and 
Doyle & Donoghue (1992). Doyle & Donoghue 
(1992) argued that this result may be an artifact 
of omitting fossils: when fossils are included in the 
analysis, this arrangement of living groups is only 
one of several that are equally parsimonious. In 
addition, several of the characters that apparently 
unite angiosperms are not unique to the group when 
fossil taxa are considered, since they also occur in 
Caytonia (flat stomata, anatropous cupules), Ben- 
nettitales (scalariform metaxylem), or both (pinnate 
sporophyll organization, integument free from nu- 
cellus). The whorled microsporophylls and tubular 
micropyle that apparently unite Gnetales are shared 
with Bennettitales. Of the characters uniting an- 
giosperms and Gnetales, opposite phyllotaxy and 
vessels are not synapomorphies if Bennettitales and 
Pentoxylon are interpolated between the two 
groups. 

Contrary to Donoghue & Doyle (19892), the 
four most parsimonious trees root angiosperms 
among paleoherbs rather than Magnoliales, with 
Nymphaeales plus monocots as the sister group of 
other angiosperms. Plesiomorphic features of Nym- 
phaeales and monocots include boat-shaped pollen 
and lack of oil cells (although the latter may not 
be valid if Acorus, which has oil cells, is basal in 
monocots, as inferred from rbcL data: Duvall et 
al, 1993). However, the conflict with the Dono- 
ghue & Doyle results is less severe than it appears: 
if angiosperms are rerooted on the branch to Mag- 
noliales, only one step longer trees are obtained 
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A one step less parsimonious tree based on the morphological data set with Magnoliales basal in 


angiosperms, showing distribution of the exine structure character. Abbreviations as in Figure 8. 


that are almost entirely consistent with trees in 
Donoghue & Doyle (1989a), except that Laurales 
are paraphyletic rather than monophyletic (Fig. 
10), and only one additional step is required to 
make Laurales a clade. The association of Mag- 
noliales with core Laurales, either as a clade or a 
basal paraphyletic group, is due to possession of 
PI type sieve-tube plastids, granular exine struc- 
ture, and a continuous tectum, all features that 
appear to be convergent when more taxa are in- 
cluded (Donoghue & Doyle, 1989a). Although 
monocots are associated with Nymphaeales rather 
than Aristolochiaceae, as in the rRNA analyses 
presented below, trees in which monocots are linked 
with Aristolochiaceae are only one step less par- 
simonious. 

This shift in rooting may be partly a function 
of the smaller sampling of taxa, but as we argued 
above in connection with the nine-angiosperm anal- 
ysis, it is also a result of treating the exine structure 
and endexine characters as unordered rather than 
implicitly ordered. Although the new morphological 
results are more consistent with the rRNA results, 
the situation does not contradict the view that the 
previous conflict between trees derived from the 
two sorts of data was a function of lower resolution 
of the morphological evidence—actually, it 
strengthens this view. Previously, morphological 
data favored a magnolialian rooting, but only weak- 
ly; now both data sets favor a paleoherb rooting, 
but the morphological data do so weakly. The fact 


that this shift followed from a rather subtle change 
in interpretation of two characters underlines the 
ambiguity of the morphological data. 

Figure 1] summarizes the bootstrap and decay 
analyses of the morphological data. The ae 
results of the bootstrap analysis are the monophy y 
of angiosperms, seen in 100% (more pem 
99.9%) of the bootstrap replicates, and, within i 
angiosperms, the association of Saururaceae 5 
Piperaceae (99%). This contradicts the tree : 
Taylor & Hickey (1992), in which weder 
and Piperaceae are distantly Vp iu 
strongest is the link between angiosperms an 4 
tales (95%). The monophyly of Gnetales is = қ 
what weaker (92%); examination of „lower: =. 
quency groupings indicates that this 18 : 
angiosperms are nested within Gnetales » i 
replicates, presumably due to features tha bs 
share with Welwitschia and Gnetum (rete E 
venation, paracytic stomata, cellular lapi 
Presumably, features of this sort are s in 
for the position of angiosperms among Gneta’ Ta 
some of the trees of Nixon et al. (1920 
Loconte & Stevenson (1990) arrangement n 
cads, Ginkgo, and conifers occurs at à ne d 
of 79%. Except for Piperales, prom Dore 
angiosperms that were seen in Donoghue 4 ’ 
(1989a) appear at much lower frequencies 
strongest being the eudicots (43%). 

The bootstrap results also bear on ! 
problem, although only indirectly. Insig 
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the регі З : 
percentage of bootstrap replicates in which each clade is found; the second number (41, etc.) indicates how many 


Steps longer trees ` 
E ees must be before some are found in which the clade no longer occurs (decays). The search of five- 


off trees was i 
AUS as ine Di FE < A 
possibility EE clades that had not decayed in that search are labeled d > 4, because we cannot rule out the 
) ey do decay in “islands” of five-off trees that were not discovered. Abbreviations as in Figure 8. 


f % 
сете frequency of angiosperm clades 
oa на ut one ог two taxa, which imply 
bine: thes not included are basal, without spec- 
groupings iF n arrangement. Ав expected, 
ils ans RAN ying that Nymphaeales and mono- 
lied 5 are most frequent, but at only 24%, 
йды.) TUS ae alone at 17%, and Mag- 
instability. of These results again illustrate the 
ues: А 1e root, while favoring the same 
The "Sai m erred from the primary analysis. 
rgt 5 numbers іп Figure 11 refer to the 
steps longer m Hat must be added (how many 
iie co | rees must be) before the group in 
date hee longer forms a clade; e.g., “42” in- 
amai А is present іп all one-off trees but 
consensus hint two-off trees (152 steps). In the 
ing within Jain trees, the only clades remain- 
and e| le are Piperales (Piperaceae 
rochodend е) and eudicots (Ranunculidae and 
"ndrales). This result reflects the unstable 


position of the root; the one-off trees include not 
only those rooted among paleoherbs and next to 
Magnoliales, but also some rooted next to eudicots, 
next to paleoherbs as a group (as in the nine- 
angiosperm analysis: Fig. 5), and next to woody 
magnoliids as a group. It also illustrates the fact 
that strict consensus trees may underestimate the 
amount of structure in the data: if a single taxon 
(or the root) “jumps” from one clade to another, 
the intervening groups collapse to a polychotomy, 
even though their other members maintain the 
same arrangement. Eudicots decay in two-off trees, 
leaving only Piperales. The arrangement of cycads, 
Ginkgo, and conifers decays in two steps, and the 
relationship between Welwitschia and Gnetum in 
four. However, angiosperms, Gnetales, the rela- 
tionship between them, and Piperales are still intact 
in the four-off trees. Five-off trees were not searched 
exhaustively because of time and memory lirnita- 
tions, but both Piperales and Gnetales decayed in 
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FIGURE 12. 






18 or more 


The two most parsimonious trees found in analysis of the rRNA data set, showing the number of 


unambiguous changes supporting each clade. Abbreviations as in Figure 8. 


several incomplete searches at this length. The 
decay of Gnetales reflects trees in which angio- 
sperms are nested within the group (cf. Nixon et 
al., 1994), since trees found by forcing angio- 
sperms together with Welwitschia and Gnetum are 
of this length (155 steps). Angiosperms and an- 
thophytes remain as clades in all five-off trees found, 
but because we cannot be certain that they do not 
decay in "islands" of trees that were not searched, 
they are labeled d > 4. This decay order closely 
parallels the relative strength of clades inferred 
from bootstrap analysis. 


RESULTS OF rRNA ANALYSES 


Analysis of the rRNA data yields two most par- 
simonious trees of 405 steps (Fig. 12), differing 
only in the relationship of Magnoliales and core 
Laurales. These trees are generally consistent with 
those derived from the whole rRNA data set (cf. 
Fig. 8), except in the exact arrangement of eudicots 
and woody magnoliids. The consistency index is 
0.58, the same as in the morphological analysis; 
the retention index is 0.66. Based on this com- 
parison, there is no reason to assume a priori that 
the rRNA data are any more or less reliable than 
the morphological data, although strong conclu- 
sions on relative consistency would be unwarranted 
because of differing amounts of missing data. 

The arrangement of cycads, conifers, and Gink- 
go differs from that derived from morphology, in 
that the group closest to anthophytes is Ginkgo 


rather than conifers. Angiosperms and Gnetales 
are united by at least 12 characters. However, it 
should be noted that conclusions on angiosperm 
outgroups are a function of the rooting of 
plants as a whole. When the whole rRNA data set 
is rooted with Equisetum and Psilotum (as in Ham- 
by & Zimmer, 1992), it is only two steps xs 
parsimonious to associate angiosperms with a " e 
consisting of Ginkgo, cycads, and conifers, rather 
than with Gnetales. The figure of 12 synapomor 
phies holds if seed plants are rooted eer 
among cycads, conifers, and Ginkgo, whic 
consider most likely. Certainly this is more сүзе, 
tent with analyses that include fossils (Crane, i е. 
Doyle & Donoghue, 1986, 1992) and wil ка 
stratigraphic record; cycads, conifers, and pee 
appear in the Late Carboniferous or poe 16 
Gnetales (or forms on the line leading to baer? 
not known before the Late Triassic (Crane, 1 
Doyle & Donoghue, 1993). 
(реак es themselves are united by at ee 
23 characters. As with previous rRNA есес 
(Hamby & Zimmer, 1992), they are rooted e 
the paleoherbs: Nymphaeales are the — ees 
of other angiosperms, followed by e 
Aristolochiaceae plus monocots. The fact с 
аге 13 unambiguous changes uniting Nymp Nym- 
and four more between the nodes exe 
phaeales and Piperales are attached may баве" 
рісіоп that this rooting is an artifact of en ae 
attraction between Nymphaeales and. агре 
groups. However, at least two considerations 
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Ficu : i 
aile x 13. Representative tree found with the rRNA data set when Magnoliales are forced to the base of 
19891) ies e paleoherbs are forced together as a clade, as in most parsimonious trees of Donoghue & Doyle 
nd Figure 10 above. This tree is 13 steps longer than the most parsimonious trees. Abbreviations as in 


Figure 8, 


et this. First, trees with the positions of Piper- 
n and Nymphaeales reversed, which are equally 
peu are only two steps longer, even 
"Weis iperales are a relatively short branch. 
ные reanalysis of the data set without Nym- 
зела results in two trees otherwise identical to 
is in Figure 12, with Piperales basal. 
CUT ehem in which alternative groups were 
iim gether with the constraints option in PAUP 
as a the view that the conflicts between 
Sees and morphological results are not seri- 
he Saga Winteraceae are interpolated 
rRNA ee and Trochodendrales in the 
M. 4 reaking up the eudicots, but forcing 
dec eid adds only two steps (if Mag- 
i € linked with core Laurales). On the other 
er vt trees rooted next to Magnoliales and 
in ene = erbs are almost equally parsimonious 
еее morphology, the magnolialian rooting 
ти SS parsimonious in terms of rRNA data. 
angios agnoliales are forced to the base of the 
Ps а the resulting trees are nine steps lon- 
MO € shortest trees (414 steps). Further- 
. » "nese trees are not closely analogous to mag- 
"rooted trees based on morphology, because 


Paleohe z 
rb groups are interpolated in various ar- 


rangements between Magnoliales and other woody 
magnoliids and eudicots. These are essentially pa- 
leoherb-rooted trees with Magnoliales alone pulled 
to the base. A more analogous tree, obtained by 
rerooting one of the most parsimonious angiosperm 
networks on the line leading to Magnoliales (Fig. 
13), is 13 steps longer than the most parsimonious 
trees. 

Experiments of this kind fail to support other 
current hypotheses on the rooting of angiosperms. 
Trees with Chloranthaceae forced to the base of 
the angiosperms (cf. Taylor & Hickey, 1992) are 
10 steps longer than the shortest trees, roughly 
the same deficit seen when Magnoliales are basal. 
Trees with Calycanthaceae basal in angiosperms 
(cf. Loconte & Stevenson, 1991) are eight steps 
longer. As when Magnoliales are forced to the base, 
paleoherbs are interpolated between the basal group 
and other woody magnoliids and eudicots. The 
shortest trees obtained by rerooting the most par- 
simonious angiosperm networks on the line leading 
to Chloranthaceae and Calycanthaceae are 11 and 
nine steps longer, respectively; the latter tree is 
relatively parsimonious because paleoherbs are a 
next-most-basal clade. Trees with eudicots as the 
sister group of other (monosulcate) angiosperms, 
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Abbreviations as in Figure 8. 


analogous to trees based on rbcL data (Les et al., 
1991; Chase et а!., 1993; Qiu et al., 1993), are 
14 steps longer than the shortest trees. 

Results of bootstrap analysis of the rRNA data 
(Fig. 14) also complement those based on mor- 
phology. Again, the monophyly of the angiosperms 
is very strongly supported (99.997%). The link 
between angiosperms and Gnetales is weaker (88%), 
but the monophyly of Gnetales is stronger (99%). 
This result supports the view that the weaker mor- 
phological support for Gnetales is due strictly to 
morphological convergences between the subgroup 
consisting of Welwitschia and Gnetum and angio- 
sperms, and it argues against trees in which an- 
giosperms are derived from (nested within) Gnetales 
(Nixon et al., 1994), 

It may be objected that molecular evidence for 
angiosperm monophyly applies only to living groups, 
leaving open the possibility that different “angio- 
sperm” lines were derived polyphyletically from 
different fossil "gymnosperm" lines. However, this 
objection is valid only if all angiosperms and their 


each 
fossil relatives are more closely related to 


other than they are to any living peo 
group, not if some angiosperm line is more 2 г > 
related to any living gymnosperm group— 80 far 
Gnetales, as assumed by most polyp ee : 
example, if some angiosperms were related to oa 
tonia and others to Gnetales, and if nl 
gave the correct relationships among жақ; ee 
angiosperms would form one branch associa ж 
Gnetales and another located one or тоге ? 
below, not a clade. 

Within angiosperms, Piperales a 
ported than they were with morpho 
they are the strongest grouping, again ver EM 
Taylor & Hickey (1992). On the other =} ak 
of the rRNA links among angiosperms dt 
flicted with the morphological results at x 
strong. The grouping of Winteraceae an gi 
odendrales, which breaks up the eudicots, hg 
at a frequency of only 54%. Although t ait 
nection between monocots and Aristolochia се 
weak (23%), monocots are linked with % 


re less well sup: 


logy (65% ). but 
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phaeales in only a negligible 1% of the bootstrap 
replicates. This may be another case where the 
rRNA data favor one alternative out of two that 
are almost equally parsimonious in terms of mor- 
phology; it requires only one extra step to associate 
monocots and Aristolochiaceae in the morpholog- 
ical analysis. The hypothesis that Aristolochiaceae 
are the sister group of monocots may be more 
plausible if Dioscoreales are basal in monocots, 
rather than alismids. 

The rooting problem can again be addressed by 
*xamining the frequency of clades containing all 
but one or two angiosperm taxa. All higher-fre- 
quency groupings of this sort imply that one or 
another combination of paleoherb groups is basal: 
Nymphaeales in 54% of the bootstrap replicates, 
Nymphaeales and Piperales in 45%, Nymphaeales 
and monocots in 25%, Piperales іп 22%, etc. In 
contrast, although Magnoliales were basal in 11% 
of the morphological replicates, they are basal in 
only 0.4% of the rRNA replicates. In other words, 
there is essentially no molecular "signal" in favor 
of the view that Magnoliales are basal angiosperms. 
This analysis also fails to support the concept that 
Calycanthaceae (Loconte & Stevenson, 1991) or 
Chloranthaceae (Taylor & Hickey, 1992) are bas- 
к Ње corresponding groupings are observed at 
Frequencies of less than 0.2%. 
M In the decay analysis (Fig. 14), Piperales and 
à © group consisting of Winteraceae, Ranunculi- 
v" and Trochodendrales are the only angiosperm 
5 es left in the consensus of one-off trees. Both 
EA ese groups decay in two-off trees. The ar- 
d gement of cycads, conifers, and Ginkgo breaks 
йры stiches steps, but angiosperms, Gnetales, 
d Wei tionship between them, and the association 
A witschia and Gnetum are still intact in five- 

trees, beyond which the analysis was abandoned. 
A aA order is generally consistent with the 
ie e strength of clades inferred from the boot- 

ар analysis, although less precisely than with the 
morphological data. 


R 
ESULTS OF COMBINED ANALYSES 


Analysis of the combined data set (Fig. 15) yields 
м ne tree of 563 steps. The consistency index 
97, which is almost identical to that in the two 
(езара, analyses (0.58); the retention index is 
d refutes one possible argument against 
morphological and molecular data, 
"amely that adding two homoplastic data sets should 
y result in more total noise. 
Xamination of Figure 15 and alternative trees 
ates graphically the complementarity of 


с 


the two data sets. In non-angiospermous groups, 
the Loconte & Stevenson (1990) arrangement is 
favored. However, the conifer-anthophyte link is 
unequivocally supported by only six characters, 
and trees with conifers and Ginkgo reversed are 
only one step longer. The strong links inferred from 
both component data sets are seen among Gnetales 
(united by at least 26 characters), among angio- 
sperms (40 characters), and between angiosperms 
and Gnetales (15 characters). 

As expected from the ambiguity of the mor- 
phological results and the relative strength of the 
rRNA results, angiosperms are rooted among the 
paleoherbs, with Nymphaeales basal. Trees with 
Magnoliales forced to the base of the angiosperms 
are 13 steps longer than the shortest tree; trees 
with Chloranthaceae and Calycanthaceae basal are 
14 and 15 steps longer, respectively. 

Despite the much greater number of rRNA char- 
acters, other results are more consistent with the 
morphological analysis, in keeping with the expec- 
tation that even a few morphological characters 
can be decisive when molecular data are ambigu- 
ous. Laurales and Magnoliales form a monophyletic 
group, as in the morphological trees, rather than 
a paraphyletic grade, as in the rRNA trees. Mono- 
cots are interpolated between Nymphaeales and 
Piperales, although it costs only one extra step to 
link them with Aristolochiaceae, as in the rRNA 
analysis. Most significantly, Ranunculidae and 
Trochodendrales form a eudicot clade— Wintera- 
ceae are dissociated from Trochodendrales and are 
instead the sister group of Laurales and Magnoli- 
ales. 

The general result of the bootstrap analysis (Fig. 
16) is that the two data sets tend to reinforce each 
other in cases where they were congruent. Angio- 
sperms (100%), Gnetales (100%), Piperales (99%), 
and Welwitschia plus Gnetum (98%) form clades 
at bootstrap frequencies similar to or higher than 
those in the separate analyses. Angiosperms and 
Gnetales are united at the 98% level, rather than 
95% in the morphological analysis and 88% in the 
rRNA analysis. This again belies the fear that con- 
flicting patterns of homoplasy will simply lower 
overall resolution, and it suggests instead that both 
analyses are detecting the same real phylogenetic 
signal. Conversely, basal seed plant relationships, 
which conflicted but were weakly supported in both 
analyses, are less resolved: conifers and antho- 
phytes are linked at a frequency of only 54%, 
rather than 79% in the morphological analysis. 

A more subtle effect is that the two data sets 
also seem to reinforce each other in one case where 
the most parsimonious trees derived from them 
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FicurE 15. Single most parsimonious tree found in analysis of the combined morphological and rRNA data set, 
showing the number of unambiguous changes supporting each clade. Abbreviations as in Figure 8. 


were not congruent. This involves the link between 
the two eudicot taxa, which is seen in 50% of the 
replicates in the combined analysis, rather than 
43% in the morphological analysis and only 12% 
in the molecular analysis. We suspect that this 
reflects the existence of “minority” rRNA char- 
acters that support the eudicots (also supported by 
rbcL data: Chase et al., 1993). Once Winteraceae 
are forced outside eudicots by the morphological 
characters, these molecular characters reinforce 
those from morphology. The existence of such ef- 
fects is a general argument for combining data sets 
(Barrett et al., 1991). 

As in the rRNA analysis, the relative support 
for clades inferred from the decay analysis (Fig. 
16) roughly parallels the bootstrap results. All the 
stronger clades decay more slowly than they did 
in either individual analysis, in keeping with the 
larger total number of characters. Eudicots appear 
to be more robust than implied by the bootstrap; 
they do not decay until four step less parsimonious 
trees. All groups retained in the four-off trees are 
still present in eight-off trees, beyond which the 
analysis was abandoned. 


CONCLUSIONS 


These exercises clearly show the utility of com- 
bining molecular and morphological data sets as 
well as analyzing them separately. This procedure 
has the potential of resolving conflicts between data 
sets even when one is much larger, presumably 


because there are minority characters in each data 
set that reflect true historical relationships (Barrett 
et al., 1991). 

The strongest results of these analyses are that 
angiosperms, Gnetales, and Piperales (Piperaceae 
plus Saururaceae, but not Chloranthaceae) e 
monophyletic groups, and that Gnetales are ES 
closest living relatives of angiosperms. We € 
that polyphyly of angiosperms can be set aside an 
other reasons examined for lack of progress 
understanding the origin of the group (cf. wee 
ghue & Doyle, 1991). In contrast, relations A 
among cycads, conifers, Ginkgo, and ante 
appear to be quite unresolved on present data, ev 
when fossil taxa are considered (Doyle & Dono- 

hue, 1992). y 
: The potentially most significant result of this 
exercise concerns the rooting of the angoa 
There appears to be essentially no rRNA m 
for the conclusion derived from morphology : 
Magnoliales are basal angiosperms, and a be: 
nolialian rooting is increasingly ambiguous in o 
of morphological data. Since there are still app 
conflicts with other molecular data sets, ìt Mes 
be premature to consider the paleoherb roo ы 
established. However, our own preliminary — 
iments in combining morphological, div 
rbcL data also give a paleoherb rooting, spec a 
between monocots and dicots. These pio ed 
not strictly comparable to the analyses P xs 
above, since so far we have only midi 
emplar species in the rbcL data set to rep 
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n DER Still, it appears that the rbcL 
i Садақ some different arrangement во 
ақы еу overwhelm the other data, even 
Sao, E more potentially informative rbcL 
Te an rRNA characters. 

ек = the rise iri among molecular data 
an EN des paleoherbs clearly deserve as 
ears cy botanists as has been paid to 
bo i A s. To discourage potential miscon- 
e р хе ould clarify what a paleoherb rooting 
ickey Ken primitive conditions (cf. Taylor & 
и dane : ). Although our trees imply that the 
the ae € would be at least semiherbaceous, 
қалы = 5 and complete lack of secondary 
аа aeales and monocots may be aut- 
NI 5 these groups. Leaves would have 
Ше с palmate venation and anomocytic 
the NIS. ку. wives Dou (1986), 
E Бағай ; etonintes of woody magnoliids would 
Sect өтесін e with Bennettitales. The fact that 
5 are near-basal might seem to support the 


View t nop 
hat small, crowded flowers with orthotropous 


ovules and no perianth are primitive (Burger, 1977; 
Taylor & Hickey, 1990, 1992). However, if taxa 
are arranged as in Figures 9, 12, or 15, it is more 
parsimonious to assume that angiosperms originally 
had flowers like those of Cabombaceae, Lactoris, 
Saruma (Aristolochiaceae), and monocots, with one 
or two cycles of three perianth parts, a trimerous 
androecium and gynoecium, and anatropous ovules. 
The syncarpous gynoecium of Nymphaeaceae and 
the tubular calyx and inferior ovary of most Aris- 
tolochiaceae would be autapomorphies, although 
the laminar placentation of Nymphaeales as a whole 
might be primitive. Stamens would be differentiated 
into filament and anther, not laminar as in woody 
magnoliids. The basal position of № ymphaeales and 
Piperales raises the intriguing possibility that the 
presence of both endosperm and perisperm in seeds 
of these orders is a primitive transitional state, not 
derived as usually assumed (D. Haig, pers. comm.). 
The conclusion that Gnetales are the closest mod- 
ern relatives of angiosperms permits a wide range 
of floral prototypes, depending on how fossil taxa 
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are arranged, from showy flowers, as in Bennet- 
titales, to simple ones, as in Gnetales. 

As stressed by Taylor & Hickey (1990), a 
(semi)herbaceous habit might help explain why the 
earliest phases of angiosperm evolution and angio- 
sperm precursors һауе been overlooked in the fossil 
record. However, this does not mean that the search 
for paleobotanical data is hopeless: paleoherblike 
precursors might be represented in the pollen, fruit, 
or seed records, and vegetative remains might be 
preserved in special facies, like Acaciaephyllum 
(a probable monocot) and Nelumbites (an aquatic 
with peltate leaves) in the Early Cretaceous (Doyle 
& Hickey, 1976). Pollen of Nymphaeales and 
Piperales is probably apomorphic in being very 
large and very small, respectively, but the two 
groups are similar in having columellar rather than 
granular structure (although the columellae are 
hard to recognize without TEM: cf. Osborn et al., 
1991) and a complete tectum. If this is the basic 
pollen type for angiosperms, it would be difficult 
but not impossible to recognize in the dispersed 
state. 

The apparent conflicts among present molecular 
data sets raise the possibility that molecular data 
are simply incapable of resolving the rooting prob- 
lem (cf. Donoghue et al., 1989), although methods 
of factoring out the effects of long branch attraction 
(as discussed by Albert et al., 1994) or discovery 
of genome rearrangements or duplications that oc- 
curred early in the angiosperm radiation (cf. Iwabe 
et al., 1989; Raubeson & Jansen, 1992) might 
permit firmer inferences. It is possible that signif- 
icant progress on the origin of angiosperms will 
require recognition of fossil forms on the long branch 
leading to the group. Potential examples include 
Phyllites (Seward, 1904), a Jurassic leaf with pa- 
leoherblike palmate venation; Triassic Crinopolles 
pollen (Cornet, 1989a), with monocotlike sculp- 
ture; and the still-enigmatic Triassic fossil San- 
miguelia (Cornet, 1986, 1989b; Doyle & Hotton, 
1991; Doyle & Donoghue, 1993). Better evidence 
on the morphology of Caytonia and glossopter- 
ids—whether or not they have anthophyte states 
in currently unknown characters, as required in 
trees where they are linked with angiosperms— 
could also have a decisive effect in choosing among 
alternative angiosperm relationships. Better recon- 
structions of primitive members or stem-relatives 
of Bennettitales and Gnetales could clarify whether 
flowers are indeed a synapomorphy of anthophytes 
or arose independently in each anthophyte line, a 
possibility raised by the existence of less flowerlike 
Late Triassic—Early Jurassic reproductive struc- 
tures related to these groups (Westersheimia, Var- 


dekloeftia, Dechellyia, Piroconites: Crane, 1988; 
van Konijnenburg-van Cittert, 1992). 

Another possibility is that evidence on the ge- 
netic control of floral development in angiosperms 
and Gnetales might indirectly distinguish among 
alternative arrangements of anthophyte groups by 
impinging upon associated scenarios of floral еуо- 
lution (Doyle, 1993). If angiosperms are basal in 
anthophytes and flowers of both angiosperms and 
Gnetales are derived from a flowerlike prototype 
(a “euanthial” scenario, as in Doyle & Donoghue, 
1986), the outer integument of Gnetales should be 
homologous with the perianth of angiosperms, and 
homologs of genes such as apetala 2 that specify 
perianth development in Arabidopsis (Coen & 
Meyerowitz, 1991) might be active during devel- 
opment of the gnetalian outer integument. On the 
other hand, trees that link angiosperms directly 
with Gnetales and place groups like Chloranthaceae 
and/or Piperales at the base of the angiosperms 
(Taylor & Hickey, 1992; Nixon et al., 1994) sug- 
gest that typical angiosperm flowers may actually 
be pseudanthia, with carpels derived from bracts 
and axillary units. If so, the outer integument of 
Gnetales might be homologous with the outer m- 
tegument of angiosperms, and its development might 
be associated with the activity of homologs of genes 
that control development of the angiosperm outer 
integument (Robinson-Beers et al., 1992). 
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APPENDIX. 


jas ORPHOLOGICAL CHARACTERS. Character definitions and 
a me follow Doyle & Donoghue (1992) for char- 
"ee 0 (except 48) and Donoghue & Doyle (19893) 
rine ang ae unless otherwise indicated. Where 
“angiosperm (fossil and d 1 2-angiosperm 
lex Ser toch ad sat pr анкер 
wee sets differ in definition of characters, defi- 
4 таң the former analysis аге given first. 
i еони (0) apical, (1) axillary. 
5 алыу buds (0) single, (1) multiple. 
dea lotaxy (0) spiral, (1) opposite or whorled, (2) 


> тей (0) ай dichotomous, (1) linear or dichoto- 
Plus cataphylls, (2) simple pinnate plus cataphylls, 


(3) pinnately compound plus cataphylls, (4) palmately 
veined (actino- or acrodromous) plus cataphylls; extant 
analysis: (0) simple pinnate, (1) linear or dichotomous, (2) 
palmately veined. States 0 and 1 of Donoghue & Doyle 
(1989a) (elliptical or obovate, secondary veins at constant 
angle or lower angle at base, vs. ovate, basal secondaries 
crowded, at higher angle) are combined under simple 
pinnate; the only taxa with Donoghue & Doyle’s state 1 
are Austrobaileya and Calycanthaceae, only one of which 
appears in each of the present data sets. 
5. Rachis (0) bifurcate, (1) simple. 
6. Laminar venation (0) open, (1) reticulate. 
7. Laminar vein orders (0) one, (1) two or more. 
8. Guard cell poles (0) raised, (1) level with aperture. 
9. Stomates (0) anomocytic, (1) mostly paracytic, (2) 
laterocytic or variable, (3) tetracytic. Core Laurales were 
scored as unknown in Donoghue & Doyle (1989a), but 
the basal state with the ingroup relationships assumed 
here (Fig. 8) is paracytic. 
10. Apical meristem (0) without tunica, (1) with tunica. 
11. Stele (0) protostele, (1) eustele with external sec- 
ondary xylem only, (2) eustele with internal secondary 
xylem. 
12. Primary xylem (0) mesarch, (1) endarch. 
13. Nodes (0) unilacunar, one-trace, (1) multilacunar 
(more than three traces from separate primary xylem 
bundles, arcuate in petiole), (2) unilacunar, two-trace, (3) 
trilacunar. The medullosan condition (many traces derived 
from one solid mass or several arcs of primary xylem, 
scattered in petiole) was treated as a separate state in 
Doyle & Donoghue (1992), but because it is uninformative 
we have eliminated the state and rescored medullosans 
as unknown; this change should have no effect on the 
results. Trochodendrales are scored as trilacunar because 
Trochodendron is polymorphic but Tetracentron is tri- 
lacunar (Cronquist, 1981). Core Laurales have various 
numbers of traces, but these are usually formed by the 
splitting of two traces (Money et al., 1950; Beck et al., 
1982), so we interpret the group as basically two-trace. 
14. Primary xylem (0) with scalariform pitting in the 
metaxylem, (1) with no scalariform pitting (coniferopsid 
type). 
Ts. Secondary xylem (0) with circular bordered pitting 
or perforations only, (1) with at least some scalariform 

iting or perforations. Scoring of angiosperms based on 
Metcalfe (1987). 

16. Vessels (0) absent, (1) present. Donoghue & Doyle 
(1989a) treated vessels in the roots only as a third state, 
but in the present data sets this occurs only in monocots. 
To preserve the unordered nature of the previous char- 
acter, monocots could be scored as unknown, but this 
would obscure the similarity between monocots and other 
groups in ability to produce the vessel cell type. To pre- 
serve this information, we һауе redefined the character 
to express this ability and scored monocots as 1. 

17. Rays (0) uniseriate ог biseriate, (1) at least some 
multiseriate. 

18. Cortical secretory structures (0) absent, (1) cav- 
ities, (2) canals. 

19. Lignin with (0) no Maule reaction, (0) Maule re- 
action (Gibbs, 1957). 

20. (0) Micro- and megasporophylls pinnately orga- 
nized, (1) microsporophylls pinnately organized, megaspo- 
rophylls simple, (2) micro- and megasporophylls simple; 
extant analysis: sporophylls (0) pinnately organized, (2) 
simple. Chloranthaceae, core Laurales, and Piperaceae 
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ТАВІЕ 1. Nine-angiosperm matrix, extant and fossil taxa. PROG = “progymnosperms” (Aneurophyton, Ar- 
chaeopteris), ELKI = Elkinsia of Serbet & Rothwell (1992) = Devonian “seed fern" of Doyle & Donoghue (1992); 
MEDU = Medullosaceae; CALL = Callistophyton; CORD = Cordaitales; CONI = Coniferales; СІМК = Ginkgoales; 
CORY = Corystospermaceae; PELT = Peltaspermum; CYCA = Cycadales; GLOS = Glossopteridales; CAYT = 
Caytonia; BENN = Bennettitales; PENT = Pentoxylon; EPHE = Ephedra; WELW = Welwitschia; GNET = 
Gnetum; MAGN = “core” Magnoliales (Magnoliaceae, Degeneria, Myristicaceae, Annonaceae); AUST = Austro- 
baileya; CHLO = Chloranthaceae; WINT = Winteraceae; EUDI = eudicots (Ranunculidae, Nelumbo, Trochoden- 
drales, Hamamelidales); ARIS = Aristolochiaceae; PIPE = Piperales (Piperaceae, Saururaceae); NYMP = Nymphaeales 
(Nymphaeaceae, Cabombaceae); MONO = monocots. ? = character state unknown; — = character not included; A 
= 0/1; B = 0/2; C = 1/2; D = 1/3; E = 2/3. 








1 2 3 4 5 6 7 8 
1234567890123456789012345678901234567890123456789012345678901234567890123456789012 
PROG @700700?7700000000?? 700700000700? 70007000070077777000007 7007? -00777771070077777777 
ELKI ??030000070000001070000701001000070007? 700070007 7777000077007? -00777771070777777277 
MEDU 10030000070070001270010701001110170007000070707?? 70000070007? -007777700700777777?? 
CALL 1003100707 1000001170010011002110071110107007 707777000007 7007? -007?777A0000777727?? 
CORD 1001700002110000012201070001211077001010700070722700000700077-00722274070027727777 
CONI 1001700000110100020201021007201000011010100000000100000000000-0072272070072727771 
СІМК 1001700000112100010201077000211000110210100000000000000000000-00227701000022772771 
CORY 7203000007201000017002011000271770111110717722777270007770077-0022227А000072777777 
PELT 2273100707222227277001001000271770110770707777222270007227277-0077770100077217777? 
CYCA ??02?00000111000120001071000711010110?10100000000000000000000 -00????010000???? m 
GLOS 1002710707127000077007000000201001111?1170??70??777000??700??-00?7???0002?7771! 
CAYT 1003710102722272222002007100201001111110717777777770007777777-00777171000077171111 
BENN 100270001? 11001012? 103010110701107 110720112220??? 70000? 2000? ?-00? 7? 701000072? ??? 7? 
PENT 1002700007272010102103010270201001110720212222222270002270077-00722701000077121777 
ЕРНЕ 1111700001112101101211170111201100210721110100010100000000000-00722701070077721774 
WELW 1111711010112101121211170111201100110721711111111100000000100-007 7? 0100002227277" 
СМЕТ 1112?110111111?1101211170711101100270?2??11111111100000000?00 00? ???1207107 ee 
MAGN 1022711111110011101002012100707001110220122100110111111100011-0001000000010040007 
AUST 1012711121112011101002017100707001110730122100110111111100001-7000001111020010000 
CHLO 101271112111201110120301171202001110730122100110711111107001-7127701111120117000 
WINT 1002711111113010101002017100707001120230122100110111111100001-000000111002001 e 
EUDI 1004211181113017101002017104707001370230122100110111111100000-010401111102 «A 
ARIS 102471110111301110100201?100707001110730122100110711111177121-?0111111100? 19 ae 
PIPE 1024711131111011101003017101707001110730122140110111111101001-1027111201420080 0 
NYMP 100471110111307070700201710070700111070122100110011111111000-70A7?1000002000 19" 
MONO 10???2111?11110?1?0100201?10070?001110?730122100110111111111020-70111101700100?? 
edie аларын ыма cies Lco лы gerer cec MEANS m 00 


x . ne 
27. Microsporophylls (0) spiral or in more ach 
whorl, (1) in a single whorl. Chloranthaceae 2% P^ 
as unknown to allow equivalence of the three-lo 


are scored as unknown to allow equivalence of their uni- 
ovulate carpels with the condition in Gnetales; however, 
Piperales as a taxon in the nine-angiosperm data set are 


scored 0, because Saururaceae have multiovulate carpels. 

21. Ovule position (0) appendicular, (1) terminal. 

22. Cupule (0) radial, lobed, (1) absent, (2) anatropous, 
(3) orthotropous unlobed; extant analysis: (0) absent, (1) 
anatropous, (2) orthotropous. 

23. Outer integument derived from two appendages 
(0) absent, (1) present. 

24. Ovules per anatropous cupule or potential homolog 
(0) several, (1) one. 

25. Microsporangia (0) terminal, marginal, or adaxial, 
(1) abaxial. Microsporangia vary from abaxial to adaxial 
in angiosperms, but we have scored them as unknown 
because of the highly modified nature of angiosperm sta- 
mens. It is problematical whether the different positions 
of the pollen sacs of angiosperms, which are unique in 
being fused lengthwise to the sporophyll, can be equated 
with conditions recognized in other seed plants. 

26. Microsporangia (0) free, (1) fused at least basally. 


droecium of Chloranthus with the whorled microspor? 
hylls of Bennettitales and Gnetales. 
d 58. Strobili (0) lacking or simple, (1) compound (и 
states in the inflorescence character of Donoghue eri 
(19892), spikes and racemes are scored aeo vah 
homologous with compound strobili, solitary ; Pt 
simple strobili; cymes are not represented, excep 
sociation with solitary flowers. 5 
"99. : Seeds (0) d (1) radiospermie, (2) рана 
ic. Omitted in the extant data set (Gnetum 15 , 
unequivocally radiospermic taxon). 
30. Integument (0) simple, (1) 
sarcotesta. 
31. Megasporangium with (0) lagenosto 
ollen chamber. С ms 
Р 32. Micropyle (0) normal, (1) tubular. Ange ұу) 
were scored as unknown in Doyle & Donoghu m 
because their ovules are so reduced, but we have 


with sclerotestà and 


me, (1) simple 
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TABLE 2. Morphological matrix, extant taxa. CYCA = Cycadales; СІМК = Ginkgo; CONI = i Я 
= Ephedra; WELW = Welwitschia; GNET = Gnetum; MAGN = “соге” Magnoliales See e dim 
Myristicaceae, Annonaceae); WINT = Winteraceae; CHLO = Chloranthaceae; CALY = росе LAUR Е 
“core” Laurales (Hortonia, Monimiaceae, Atherospermataceae, Siparunaceae, Gomortega, Hernandinceso Laura- 
сеге); SAUR = Saururaceae; PIPE = Piperaceae; ARIS = Aristolochiaceae; NYMP = Nymphaeales [Ny ipliaen сане 
Cabombaceae); RANU = Ranunculidae; TROC = Trochodendrales; MONO = monocots. ? = character state Блейк 
= = character not included; A = 0/1; B = 0/2; C = 1/2; р = 1/3; E = 2/3. 


Oo COMU 
T ——M—M—M————————————————————— 


1 2 3 4 5 6 7 8 

бл САҒА АЕ 3507 8401774567890120 RONL асет мМ 496040121967020 1A 356759012 
dui -?00-0-000- -10001201000-1000-1-0-00- - -00-00000000000000000000000?? ? ?010000? ??? ???? 
E -001-0-000- 21000100000 -1000-1-0-00-- -00-00000000000000000000000??? ?010000? ??????1 
Rar eea O -010002000001007 -0-0-00 ---DO-00100000 100000000000007 ARDOR 7777771 
Ki Ree pr A-T E икен og dae Ed 
ЕТ -111-1-010--21011210101-0111-0-1-00---11-11111111100000000100000? ?? ?010000?? ?? ???? 

-110-1-011--11?11010101-0?11-0-1-01---1?-111111111000000007000007???120?10???????1 


ня -020-1-111--00111011010-7100-0-0-10---10-221001101111111000117000100000001000000? 
ENT -009-1-111--30101011010-7100-0-0-10---20-22100110111111100001700000011100200111001 
eee 21--20111017020-7171-0-0-10---20-22100110711111107001771277011111201170000 
ve :010-1-111--20111011010-7100-0-0-17---20-22100170211111700111000000010070210300072 
ДІҢ -018-1-111--20111017010-7100-0-0-11---10-2210011071111110901110700011107127127000? 
ae 7022-1-131--10717011020-7101-0-0-10---20-22100110711111101001010271112010200010111 
Mcr 1--1011101?020-?101-0-0-10- - -?0-22110110111117102?017102?1112011707210111 
ms -922-1-101--30111011010-7100-0-0-10---20-221001107 111111771217 701111111007 10010000 
er 1-101- -30707071010-7100-0-0-10---CO-22100110011111111000070477100000200011011 
codec ana 01--00111011010-2100-0-0-12---20-22100110111111107000101017211717200771000 
DE -121--30101011010-?101-0-0-12---20-2210011021111110000007120011111020001070? 

?2-1-1?1--10?1?011010-?10? -0-0-10---20-22100110?11111111020070111101?00100???00? 


с ЗВЕР ers 


39. Exine structure (0) spongy alveolar, (1) honey- 


them as normal since the i i 
. y certainly have nothing closel 
E pon to the siete аанак: мы. 4 
4, eas (0) not vascularized, (1) vascularized. 

reel) i lar cuticle (0) thin, (1) thick; extant analysis: 
Cle онй integument (0) fused most of the way 
eg se us and (1) free nearly to the base, which 
bees ave the same distribution but is better docu- 

35. р, cn taxa. 
ais (3) a with (0) tetrad scar, (1) sulcus, (2) no ap- 
n (0 ree colpi; extant analysis: (0) sulcus, (1) no 
tetrad 5A 2) three colpi. Conifers are scored as having a 
minit E i the nine-angiosperm data set, based on 

ape The à representatives, but as sulcate in the extant 
+ eds sc isulculate pollen of Calycanthaceae is scored 
арша, s allow derivation from either sulcate or 

Еторе they were scored as state 1 in Donoghue & 
ё а), but this state included sulcate, inapertur- 

wp sulculate. 
Doyle Қа ауа (0) radial, (1) bilateral. As іп 

fia its десі ue (1992), we score Gnetum as unknown, 
ipn ie len has global rather than radial or bilateral 
E te е also score angiosperms with radial pollen 
to be eri not only because their symmetry is unlikely 
Pollen of gous with radial symmetry in the sporelike 
correlated peat, seed plants, but also because it is 

24 it „other characters included in the matrix: 
Uninfo pe in Winteraceae, three colpi in eudicots. 
е іп the extant data set. 
ollen (0) nonsaccate or subsaccate, (1) saccate. 

ceus structure (0) eusaccate (alveolae detached 


from 4 

пе 

lectum лац (1) protosaccate (alveolae continuous from 
to nexine), 


comb alveolar, (2) granular, (3) columellar; extant anal- 
ysis: (0) alveolar, (1) granular, (2) columellar. Unordered 
for reasons discussed in text. Nymphaeales are scored as 
uncertain because Osborn et al. (1991) reported that 
Cabombaceae are columellar. 

40. Exine striations (0) absent, (1) present. Doyle & 
Donoghue (1992) scored Gnetum as unknown, on the 
grounds that its tectum is so reduced that any striations 
would have been lost; this is now supported by ultrastruc- 
tural observations that the tectal spines of Gnetum re- 
semble striations of Ephedra and Welwitschia (Gillespie 
& Nowicke, 1992). 

41. Megaspore tetrad (0) tetrahedral, (1) linear. 

42. Megaspore wall (0) thick, (1) thin, (2) absent (or- 
dered). Ordering of this character is in keeping with the 
treatment of thin and absent as one state in Doyle & 
Donoghue (1986, 1992). 

43. Microgametophyte with (0) more than four nuclei, 
(1) four nuclei, (2) three nuclei. In Doyle & Donoghue 
(1986, 1992), the three-nucleate state of angiosperms 
was omitted and angiosperms were scored as unknown 
(to allow derivation from either state), but this state is a 
potential synapomorphy of angiosperms in the present 
data sets. 

44. Sperm transfer (0) zooidogamous, (1) siphonoga- 
mous. Conifers are scored as zooidogamous in the nine- 
angiosperm data set, based on lack of a sulcus in primitive 
fossil representatives (cf. Doyle & Donoghue, 1992), but 
siphonogamous in the extant data set. 

45. Megagametophyte (0) monosporic, (1) tetrasporic. 

46. Egg (0) cellular, (1) free-nuclear. 
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47. Early embryogenesis (0) free-nuclear, (1) cellular. 

48. Fertilization (0) single, (1) double. See discussion 
in text. Calycanthaceae lack double fertilization (Loconte 
& Stevenson, 1991), but because this is associated with 
apomixis (Davis, 1966) we score them as unknown, 

49. Embryo (0) without feeder, (1) with feeder. 

50. Seed germination (0) hypogeal, (1) epigeal. An- 
giosperm data from de Vogel (1980) and Endress (1983); 
contrary to Loconte & Stevenson (1991), Endress reports 
that germination of Austrobaileya is epigeal. 

51. Companion cells in phloem (0) absent, (1) present. 

52. Microsporangia (0) various, (1) in two pairs. 

53. Endothecium (0) absent, (1) present. 

54. Pollen germination (0) in pollen chamber, (1) on 
stigma. 

55. Megagametophyte (0) large, (1) eight-nucleate. 
The larger tetrasporic embryo sacs of Piperaceae are 
conservatively scored as unknown, although they show 
similarities with the eight-nucleate type. 

56. Endosperm (0) absent, (1) present. 

57. Radicle (0) persistent, (1) replaced by adventitious 
roots. 

58. Habit (0) woody, (1) herbaceous. Groups with in- 
terfascicular cambium not producing normal secondary 
xylem are scored as unknown. 

59. End-wall pits or vessel perforations (0) multiple, 
(1) simple. 

60. Sieve-tube plastids (0) starch, (1) PI type, (2) PII 
type. See text for addition of PII type, scoring of Mag- 
noliales (Behnke, 1988). Core Laurales were scored as 
unknown in Donoghue & Doyle (1989a), but the basic 
state with the ingroup arrangement assumed here is PI. 

61. Oil cells (0) absent, (1) present. 

62. Benzylisoquinoline alkaloids (0) absent, (1) pres- 
ent. Uninformative in the nine-angiosperm data set. 

63. Stipules (0) absent, (1) adnate-axillary. Other stip- 
ule types are scored as unknown. 

64. Chloranthoid teeth on leaf margins (0) absent, (1) 
present. 

65. Perianth (0) more than two whorls (or spiral or 
chaotic), (1) two whorls, (2) absent. Donoghue & Doyle 
(1989a) scored Trochodendrales as having one whorl, but 
because this state is autapomorphic in the present data 
set we have rescored the group as unknown. 

66. Perianth symmetry (0) various, (1) at least calyx 
trimerous. 

67. Stamen number (0) various, (1) multiples of three. 

68. Stamens (0) laminar, (1) with well-differentiated 
filament. 

69. Pollen (0) boat-shaped, (1) globose. Groups with 
saccate and sporelike pollen scored as unknown because 
it is unclear whether their shape conditions can be com- 
pared with those in nonsaccate, basically monosulcate 
groups. 

70. Pollen size (0) large (> 50 um), (1) medium, (2) 
small (« 20 um) (ordered). Previously this character was 
used only in the angiosperm study (Donoghue & Doyle, 
1989a). We have scored taxa with saccate pollen (except 
Caytonia) as uncertain (0/1), for two reasons. First, 
saccate pollen tends to be conspicuously larger and more 
massive than pollen of nonsaccate groups, suggesting that 
there is a functional correlation between presence or ab- 
sence of sacs and size. If so, scoring saccate groups on 
the angiosperm-based size scale of Donoghue & Doyle 
might excessively weight transitions between saccate and 
nonsaccate and produce spurious groupings among non- 


angiospermous taxa. Second, it is possible that the more 
appropriate comparison is between size of nonsaccate pol- 
len grains and size of the central body (nexine) in saccate 
pollen, which is often much smaller than total grain size, 

71. Tectum (0) continuous or finely perforate, (1) fove- 
olate-reticulate. We have modified the definition of state 
1 from semitectate-reticulate in Donoghue & Doyle (19895) 
and rescored Austrobaileya as having this state, since its 
relatively large foveolae and rounded muri seem more 
comparable to the sculpture of Chloranthaceae and other 
reticulate groups than the much more finely perforate 
sculpture of groups such as Magnoliales. 

72. Aperture membrane (0) (nearly) smooth, (1) sculp- 
tured. Conifers are scored as unknown in the nine-angio- 
sperm data set because a sulcus is absent in primitive 
forms. Scoring of Aristolochiaceae is based on SEM photos 
of Saruma pollen kindly provided by Long Huo (Guang: 
zhou); scoring of Piperales is based on Bornstein (1989) 
and Doyle & Hotton (1991). 

73. Supratectal spinules (0) absent, (1) present. Scor- 
ing of Piperales is based on Bornstein (1989) and Doyle 
& Hotton (1991). 

74. Endexine (0) thick, laminated, (1) absent, (2) thin, 
nonlaminated, except under apertures. Unordered for rea- 
sons discussed in the text. Following Chlonova & Surova 
(1988), we have rescored Chloranthaceae as having end- 
exine, not unknown as in Donoghue & Doyle (1989a). 

75. Hypanthium (0) absent, (1) present. же.” 

76. Оушев per carpel (0) several, (1) опе apical. / 
third state in Donoghue & Doyle (1989a), one —— 
represented only in Piperaceae, which are therefore scor 
as unknown. . 

77. Fruit (0) dehiscent, (1) berry, (2) drupe (with e 
docarp), (3) dry indehiscent. As noted by Loconte 
Stevenson (1991), the spongy, indehiscent or ie Н 
dehiscent fruits of Nymphaeaceae (Cronquist, 1981), 
scored as dehiscent in Donoghue & Doyle (19894), are 
better characterized as berries. Winteraceae were Wes 
as unknown in Donoghue & Doyle (19893), for bem 
dehiscent or berries, but because the dehiscent conditi 
is restricted to Takhtajania, which appears to inb 
within the family (Vink, 1988), we have rescore 
as having berries. bae 

78. Testa (0) multiplicative, (1) non-multiplicative. 

79. Exotesta (0) e aree 

80. Tegmen (0) normal, sclerotic. 

81. Nutritive tissue (0) ule ae 

lus perisperm. Megagametophyte tissue : 
sidered a ыға vem jam this xm be redundant wi 
endosperm formation (character 56). 2 12- 

82. Chromosome number (0) n = 6-8, (1) i o i 
16. Donoghue & Doyle (1989a) defined state жы» 
12-19, but because numbers above л = 12 ure 
certain and potentially complex origin, we ках кет: 
the states and scored taxa with n > 16 as 2 10-11 
third state in Donoghue & Doyle (19892), ” here 
i which are t 
is represented only by Calycanthaceae, 
fore scored as unknown. Data оп non-angiospermous 
from Ehrendorfer (1976). 

RiBosoMAL RNA CHARACTERS. 1-86: 
the 18S subunit, in terms of positions іп soy. 

1: 93; 2: 120; 3: 131; 4: nucleotide 189 10: 190: 
soy; 5: 132; 6: 134; 7: 176; 8: 181; 9: 191-1925 

11: nucleotide not present іп soy eres 17: 240 
12: 200; 13: 220; 14: 222; 15: 236; 16: 29% 

18: 241; 19: 242; 20: 244; 
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TaBLE 3. Ribosomal RNA matrix. Order of taxa as in morphological matrix. R = A/G; Y = C/T; M = A/C; 
K = G/T; S = C/G; W = A/T; Н = A/C/T; B = С/С/Т; V = A/C/G; D = A/G/T; М = A/C/G/T. 
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ACCTCGCTTTTSCGCGTBDTCGWAYTGCTMCASGA TAGCYGACACCATATTGCCCRGTTTTTGCCTGTGCGTAATACATTAGCTTAG 
ACCCTACTTC2CCGCCTGCTCGTACCGCTGCACGATA? 2? TGATACCGT? TTGCTCGTTCGTCGCCT? 2222 2AATAAATTAGCTTA? 
GYCYWRCTCTTTCGCCTWCCCGWACCRYTTTASGATAGCCGACACCCTGGTGCCCGTTCCTCGCTT? 2 GCCGAATACRTTAGYTTAG 
GCCCTAGTTC?TCGCCTTGTCGAATCGTCACWYKRCHGTCGATTCCCCGTCG2CTCATTCCGGATTGCGTGGTACGGGCCAACTTGG 
GCCCTWGTTC?TCATCGTGTTGAGCCGTCACACTGCA 2AGGGTTTCCCGGTGTCTCGCAGTTGATT? 22 22? TACGCGCCAATTTGG 
?CCCTAGTCT?TCGTTTTTTCGAGG?ACTTCAGTGCA 2 TCGACTTCCCGGTGTCTGACGCTCGACT? 22 222 2A ? GGGTCAATTTGG 
AGACYGACTCGTYGCTCGCCTATATCGCTCYWKGACTGCCGATTCCATATCGCYYGTGACCCACTCTC?2 GGTAGTGCGTTGGCCCAT 
AGT?CAACTC? TTACTTGCCTAAACCGCTC2WYKRCWGACAATTCTATATCGYTTGTCACCCACTCTCTKRTAGTGCGTTGGCCCAY 
AGT?CGACTC? TTGCT?2GCCTATACCGCTTCTCGACT?TCGATTCCATATCGCTTGT2ACCCACTC?2? ? ? TAGTGCGTTGGCCTA? 
?GTCT?CCTC?TTGCTCGCCTATACCGTCTCTTGACTGACGATTCCATATCGCTTGCGACCCACTC? ? T2 GTAGTGCGTTGGCCCAT 
AGT?CTACTC?TTGCTCGCCCATAYYGCTCCGTGACTTACGATTCCATATCGCTTGTGAYCCACTS ? 2? TGTAGTGCGTTGGCCCA ? 
AGC?TGCCTC?TTACTTGATT?TAC?GTCT? TCGACATACGGTTCCATATCTCCCGATCCCCACTCTCGTTTAGTGCGTTGGCCCAT 
AGT?TRCCTCGTTAYTYKMMYATATCGTTTTKCKACW2ACGATTCCATATCTCCTGYCWCCCACTY ? TGTTTAGTGCGTTGGCYYAC 
AGC?TGMCTC?YTACTTGYYTATACTGYTTCGCRACT2ACGATTCCATAKCGCTTGTGACCCACTCTCYGATAGTGCGTTGGCCCAT 
AGC?TRCCTC? CCRCATGTTTGTACCGYTYCWCAACTGATGAYTCCATATCGCCTGWMCCCCRCTC? ? GGGAAGTGMGTTGGCTTAK 
AGT2AGACTC?TTATTTGCCTATA?2222??2 TTGACTGCCGATTCCMTATCGGTTGTACCCCACTCTCTGATAGTGAGT TGACCTAT 
AGT?CGACTC?TTGCTTGCCCATACCAC?? TATGACT? CCAATTCTATATCGCCTGTAACCTACTCTCTGATAGTGCGTTGGCCCAT 
AKCKSGMCTC?YYGCNTGMYBAYAYHGYTCCGCGACTGACGAYWCCATAYCGCTCGYRWCCCGCTCTCKGGTARTRCGTTGGYCCAG 


3 0 1 2 3 4 5 6 7 
890123456789012345678901234567890123456789012345678901234567890123456789012345678901234 
TCGCTCTATATTAG?CCCCRATCGCCTATTWACCTCGCGYCTATGY CTGGCRYTCTWYRTWCGYSSGGGTTAGGGTGATGGGAGAGR 
222??? TATACCAG? CCGCAATCGCCTATTAACCCCGTGTTTATGCCCGGCGCTCTACACCCGTCGCGGTTAGGGTGATAGG?G7AG 
TCGTYCCATACYARCCCGCAATCGCCTRTTVACYYCGTGTCTATGCCTGGYGCYCYHCATCCGYTCGGGTTAGKGTGATAGGGGARG 
?TCTTCCATGCCAACCCGCGACTGCTTATTAAC? TCGCGCCTATRCTCGCCTGCTAACGCCCGCCTCGGTTAGAGCGTCTGGGGAGA 
ATCCAGTGTATTCACCCGCAACCACTTATAAATCACGTGCCTATGCCCACTTGCCTCCGTCAGSTTGGG? T TGAACGTCTGA 7A GGG 


TIGTAATGCGCC2ATCCTCAGTCGTCGGACKTCGCTGCGC2 22222222 222222222222222 2222222??? 2??? ? ?A6G6AGG 
TIGTATTACGCCAATCCTCAGCTGTCGGWCGGMGCTGCGCTCGCA?CCGCCGCATTCCGCGCGCMTCSGCCGRTGCRATAAGGGAGG 
TTGAATTACGCMAA? 2TTCGGTCGTCGGACGGCGCTGCGC2222 2 2 CCGCCGTATT? CGCGCGTYYCGGCCGGTGTRATAAGGGAGG 
TCGTAGTGCGTCAR?TYTCGGTTGTCGGACGRCSYYAYAYWCGCATTCGCCGYATTY CRCGCAYTYVTSCCGATGCRATAAGGGAGG 
TCGATTTACGCCAATCTTCGGTCGTCGGACGGCGCTGCGCACGCA?TCATTGCATTCCGTGTGCTTCGGCCGGTGTGATAAGGGAGG 
TCGAAATGC?CAAGGCTTTGGGTGTCGAACGACGCTGCGC? 222222222222??? ? 2 GCGCGTT? CGGCCGGTGTGATAAGGGAGG 
TYGHNGTGCRCCWAYCCTCGGYCGYCGGACGGCGCYGCGCHCGSG2 CCGCYGYATTYYGCGCGTYKCGGCCGGTGY GATARGRRARG 
MISES и Oe EON CRE GNC CULES uot cans DUAE ee ae 


oral 245; 22: 247; 23: 251; 24: 257; 25: 263; 26: 91: 759; 92: 769; 93: 784; 94: 790; 95: 791; 96: 
215; 27: 276; 28: 280; 29: 281; 30: 282; 797; 97: 830; 98: 834; 99: 865; 100: 866; 
31: 287; 32: 339; 33: 347; 34: 353; 35: 368; 36: 101: 904; 102: 1602; 103: 1612; 104: 1620; 105: 
m 37: 496; 38: 542; 39: 936; 40: 1042; 1621; 106: 1624; 107: 1637; 108: 1639; 109: 1650; 
l: 1050; 42: 1055; 43: 1063; 44: 1065; 45: 1075; 110: 1651; 
46: 1076; 47: 1085; 48: 1086; 49: 1096; 50: 1100; 111: 1656; 112: 1662; 113: 1663; 114: 1683; 115; 


_ 51: 1241; 52: 1245: 53: 1300; 54: 1355; 55: 1357; 1705; 116: 1708; 117: 1712; 118: 1716; 119: 1723; 
56: 1358; 57: 1363; 58: 1364; 59: 1365; 60: 1366; 120: 1731; А. Рах 
61: 1370; 62: 1371; 63: 1372; 64: 1376; 65: 1404; 121: 1757; 122: 1758; 123: 1760; 124: 1764; 125: 
6: 1411; 67: 1503; 68: 1514; 69: 1526; 70: 1527; 1777; 126: 1783; 127: 1796; 128: 1949; 129: 1950; 
,,11:1528; 72: 1534; 73: 1555; 74: 1564; 75: 1566; 130: 1951; 4 
16: 1568; 77: 1573; 78: 1606; 79: 1613; 80: 1666; 131: 1958; 132: 1960; 133: 1961; 134: 1962; 135: 


в6.1:1668; 82: 1677; 83: 1724; 84: 1729; 85: 1735: 1969; 136: 1971; 137: 1973; 138: 1980; 139: 1982; 
: 1747. 140: 1983; - 
87-167: Characters from the 26S subunit, in terms 141: 1986; 142: 1991; 143; 2000; 144: 2001; 145: 


of Positions in rice. 2028: 146: 2031; 147: 2032; 148: 2037; 149: 2040; 
87: 740; 88: 741; 89: 748; 90: 750; 150: 2041; 
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151: 2056; 152: 2057; 153: 2058; 154: 2059; 155: 
2060; 156: 2061; 157: 2066; 158: 2072; 159: 2073; 
160: 2077; 

161: 2078; 162: 2085; 163: 2086; 164: 2088; 165: 
2089; 166: 2098; 167: 2102. 

168-174: Insertion-deletion events, where gap = А. 

168: indel 1 between 131-132 in soy; 169: indel 2 
at soy 454; 170: indel 8 between 1526-1527 in soy; 

171: indel 9 between 1593-1594 in soy; 172: indel 
10 between 768-769 in rice; 173: indel 11 between 
1746-1750 in rice; 174; indel 12 between 1756-1758 
in rice. 


ADDENDUM 


VOUCHERS OF SPECIMENS USED IN rRNA ANALYSES 


Aristolochia sp. (probably A. gigantea Mart. & Zucc.): 
Towa State Univ. greenhouse, provided by J. Wen- 
dell, Z-94-1, 1/5, 

Arundinaria gigantea (Walter) Muehlenb.: Baton Rouge 
Parish, coll. L. Sims, Z-12-86, LSU. 

Asimina triloba (L.) Dunal: Burden Plantation, Baton 
Rouge, coll. M. Bowen & C. Knaak, Z-16-89, LSU. 

Barclaya longifolia Wallich: Suwanee Labs, Lake City, 
FL, provided by D. Bryne, Z-8-89, LSU. 

Brasenia schreberi J. Е. Gmel.: Golden Ranch Farm, 
Gheenes, LA, provided by B. Crain, Z-9-89, LSU. 

Calycanthus occidentalis Hook. & Arn.: Hilltop Arbo- 
retum, Baton Rouge, coll. C. Knaak, Z-7-89, LSU. 

Chloranthus spicatus (Thunb.) Makino: Davis Botany 
greenhouse B81-804, Doyle 94-2-09.3, DAV. 

Colocasia esculenta (L.) Schott var. antiquorum (Schott) 
Hubb. & Roeder: LSU campus, coll. L. Sims, Z-1- 
85, LSU. 

Cycas revoluta (Thunb.); LSU campus, Z-2-85, LSU. 

Echinodorus corditollus Griseb.: coll. F. Givens, Z-2-88, 


Encephalartos ferox Bertol. f.: Univ. Illinois greenhouse, 
provided by D. Nickrent, Z-94-2, US. 

Ephedra distachya L.: Davis Arboretum A65-888, Doyle 
94-2-09-2, DAV. 

Ephedra tweediana C. A. Mey.: Davis Arboretum A67- 
620, Doyle 94-2-09-3, DAV. 

Euryale ferox Salisb.: D. Les s.n., CONN. 

Ginkgo biloba L.; LSU campus, coll. L. Sims, Z-3-85, 
LSU. 

Gnetum montanum Markgr.: Davis Botany greenhouse 
B70-116, Doyle 94-2-09-1, DAV. 

Бауды sp.: Mt. Kogi, New Caledonia, L. Thien 600, 


Hosta japonica Tratt.: Baton Rouge, coll. L. Sims, Z-17- 
87, LSU. 

Juniperus ashei Buchholz: Travis Co., TX, coll. J. Drost, 
2-22-87, LSU. 

Liriodendron tulipifera L.: Baton Rouge, coll. C. Knaak, 
Z-5-89, LSU. 

Magnolia grandiflora L.: LSU campus, Z-4-85, LSU. 


Najas guadaliensis (Spreng.) Magnus: coll. C. Knaak, 
Z-1-89, LSU. 


Nuphar luteum (L.) Sibth. & Sm. subsp. macrophyllum 
(Small) Beal: Jefferson Co., MO, P. H. Raven 27204, 
MO. 

Nymphaea odorata Ait.: Tammany Parish, LA, coll. F. 
Givens, Z-21-87, LSU. 

Persea barbonia (L.) Sprengel: Hilltop Arboretum, Baton 
Rouge, coll. C. Knaak, Z-6-89, LSU. 

Piper nigrum L.: Iowa State Univ. greenhouse, provided 
by J. Wendell, Z-94-3, US. 

Potamogeton sp.: P. Hoch 3439, MO, 

Ranunculus acris L.: coll. L. Sims, 2-3-86, LSU. 

Sabal minor (Jacq.) Pers.: Burden Plantation, Baton 
Rouge, coll. J. Drost, Z-20-87, LSU. 

Sagittaria lancifolia L.: Baton Rouge, coll. E. Jupe & 
E. A. Zimmer, Z-13-86, LSU. 

Saruma henryi Oliver: US National Arboretum 49482, 
coll. J. Kress, Z-94-4, US. 

Saururus cernuus L.: San Gabriel Parish, LA, coll. R. 
Chapman, Z-3-88, LSU. 

Trochodendron aralioides Sieb. & Zucc.: Taiwan, S.-M. 
Chaw 189, HAST. 

Zamia floridana A. DC.: Univ. Illinois greenhouse, pro- 
vided by D. Nickrent, Z-94-5, US. А 

Zamia ottonis Miq.: Univ. Illinois greenhouse, provided 
by D. Nickrent, Z-94-6, US. 


SPECIMENS WITHOUT KNOWN VOUCHERS 


Avena sativa L.: seeds provided by S. Roux, Univ. Texas. 
Cabomba caroliniana A. Gray: San Marcos River, TX, 
rovided hy E. Schneider. 
CRIME japonica (L.f.) D. Don: Mellingberg дн 
Drimys winteri Forster & Forster f.: Berkeley Вогалі 
Garden 45.307, provided by J. Affolter. 
Hordeum vulgare L.: cultivar "Himalaya," seeds pro- 
vided by M. Saghai-Maroof. Т ач 
Oryza sativa L.: cultivar Lamont, Louisiana Rice 
search Station. eec 
Peperomia sp.: Univ. Illinois greenhouse, provided by 9. 
Nickrent. 2; 
Pinus taeda L.: seeds provided by O. Stubbs, Louisiana 
Dept. Wildlife & Forestry. ema 
Pistia stratoides L.: provided by P. Hoch, Missouri 
tanical Garden. : EU 
Saccharum officinarum L.: line Cavengerie, pro 
K. Damann, LSU. 
Sorghum bicolor (L.) Moench: cultivar Tr P A 
provided by M. Thomas-Compton, Univ. as 
Tasmannia lanceolata (Poir.) A. C. Sm.: pri : 
tanical Garden 60.0052, provided by J. char 
Tripsacum dactyloides (L.) L.: line Tp 1 12 of i 
provided by K. Newton, Univ. Missouri. imn 
Triticum aestivum L.: line HW3022, Rohm an 
seeds, provided by S. G. Bartlett, LSU. дәлдігі 
Welwitschia mirabilis Hook. f.: Huntington en 
Gardens, J. Folsom 13, provided by E. pin 
Zea mays L.: cultivar B73, Pioneer HiBred " 


